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Computer  simulation  of  transient  heat  transfer  from  coated  superconductors 
to  liquid  helium  have  revealed  that  coating  a  superconductor  with  certain 
new  J”® rial s ,  called  Laketites,  instead  of  traditional  insulation,  such 
as  GE7031,  substantially  improves  the  ability  of  the  superconductor  to 
withstand  transient  heat  pulses  without  making  the  transition  to  the 
normal''  i.e.  non-superconducting  state.  In  the  best  cases,  the  length 
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of  time  that  constant  heating  can  be  withstood  is  improved  twenty  fold. 
This  report  shows  that  the  most  important  property  of  these  new  materials 
is  their  increased  thermal  conductivity  relative  to  current  materials. 

The  highest  thermal  conductivity  materials  designated  SC-2,  SC-3  are  the 
most  desirable  for  further  research  and  development.  The  increased 
specific  heat  of  the  Laketites  is  desirable,  but  not  crucial.  Further¬ 
more,  the  thickness  of  the  coating  has  little  effect  on  the.  improved 
ability  of  the  superconductor  to  withstand  heat  pulses. 

This  report  recommends  a  vigorous  program  of  further  development 
of  the  Laketites.  A  complete  theoretical  background  of  the  computer 
simulation  and  examples  of  the  programs  used  are  also  included  in 
this  report.  „ 
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I  Summary 

The  objective  of  this  investigation  was  the  computer  simulation  of  tran¬ 
sient  heat  transfer  from  superconductors  coated  with  the  newly  developed  Lake- 
tite  coatings  to  liquid  helium  in  order  to  test  the  ability  of  these  coatings 
to  improve  the  heat  transfer  process  at  low  temperatures  and  to  make  recom¬ 
mendations  as  to  which,  if  any,  of  the  three  classes  of  materials  should  be 
selected  for  further  research  and  development.  A  secondary  objective  is  the 
investigation  of  the  changes  in  heat  transfer  characteristics  due  to  variation 
of  the  thickness  of  the  coatings. 

A  computer  program  using  the  SCEPTRE  code  was  developed  to  perform  the 
computer  simulation.  An  example  of  an  actual  SCEPTRE  program  is  given  in 
Appendix  I  along  with  explanatory  notes.  Since  each  material  has  different 
thermal  properties  the  program  must  be  changed  in  specified  places  to  reflect 
these  changes  in  properties.  Critical  to  the  operation  of  the  program  was  a 
Fortran  subroutine  which  is  given  in  well -documented  form  in  Appendix  II.  A 
computer  simulation  is  only  accurate  to  the  extent  that  the  equations  that  are 
used  accurately  reflect  the  physical  reality  of  the  process  that  is  occurring. 

The  heat  transfer  theory  that  underlies  the  equations  used  in  the  computer 
simulation  and  all  the  assumptions  and  approximations  used  to  derive  these 
equations  are  presented  in  detail  in  the  theory  section  of  this  report.  All 
of  the  equations,  approximations,  and  assumptions  are  referenced  to  existing 
heat  transfer  literature  to  provide  a  justification  for  their  use. 

The  computer  program  thus  developed  was  used  to  test  SC-1C,  SC-2B,  and 
SC-3A  Laketite  coatings  with  bare  copper  and  GE7031--a  typical  present-day 
insulation— used  as  comparison.  The  superconductor  was  represented  by  a  5  mil 
thick  flat  copper  plate  to  imitate  the  copper  outer  sheath  usually  present  on 
superconducting  wire.  Constant  heating  with  power  levels  of  .05,  .1,  .2,  .5, 

1,  2,  3,  4,  5,  10,  20,  and  50  W/cm^  was  used  to  span  the  range  of  heating 
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power  likely  to  be  encountered.  Constant  heating  is  produced  by  the  eddy 
current  process  due  to  changing  magnetic  fields.  Coating  thicknesses  of  1, 

2,  4,  10,  and  20  thousandths  of  an  inch  (mils)  were  used.  The  time  of  heating 
was  varied  from  .001  seconds  to  .1  seconds.  Despite  these  particular  choices, 
the  program  was  designed  so  that  it  could  accept  any  material  whose  thermal 
properties  were  well  characterized*,  any  shape  of  heating  pulse,  length  of 
heating,  and  level  of  heating,  as  well  as  any  changes  in  bath  temperatures. 

From  the  voluminous  results  of  the  computer  simulations,  three  principle 
conclusions  emerged.  1)  The  Laketites  provide  a  substantial  improvement  in 
the  thermal  response  characteristics  over  present  materials.  Thus,  a  vigorous 
program  of  further  research  and  development  on  these  materials  should  be 
instituted.  2)  The  most  important  parameter  for  a  coating  to  have  is  a  high 
thermal  conductivity,  with  specific  heat  being  less  important.  Thus,  the  SC-2 
and  SC-3  materials  appear  to  be  the  most  promising.  3)  The  thickness  of  the 
coating  is  not  a  vital  parameter  and  should  be  governed  by  mechanical  consid¬ 
erations. 

Recommendations  for  further  research  and  development  of  these  coatings 
and  developments  of  better  computer  simulations  are  made. 


II  Introduction 


The  Air  Force  has  long  been  interested  in  applying  the  phenomenon  of 
superconductivity  for  the  purpose  of  lightweight  production  and  storage  of 
energy.^  Superconductors  have  the  disadvantage  that  they  must  be  refrigerated 
below  a  transition  temperature  in  order  to  operate.  If  a  section  of  super¬ 
conductor  heats  up  above  its  transition  temperature,  that  section  becomes  a 
"normal"  high  resistance  section  of  wire.  The  continuing  passage  of  electrical 
current  through  the  resistive  section  produces  additional  ohmic  heating  which 
raises  more  of  the  superconductor  above  its  transition  temperature,  leading 
to  a  runaway  reaction  that  disables  the  superconducting  device.  This  prob¬ 
lem  is  particularly  severe  for  devices  that  are  pulsed,  since  pulsed  operation 
creates  rapidly  changing  magnetic  fields  that  produce  large  eddy  current 
heating  in  the  copper  or  bronze  matrix  that  surrounds  the  superconducting 
portion  of  the  wire.  The  shorter  the  pulse,  the  worse  this  problem  becomes. 

As  superconducting  energy  storage  devices  may  well  be  operated  in  the  pulsed 
mode,  any  way  to  increase  the  ability  to  remove  or  absorb  the  heat  produced 
during  a  pulse,  thus  preventing  the  superconductor  from  becoming  "normal" 
is  desirable. 

In  1980  Lake  Shore  Cryotronics,  a  concern  that  manufactures  cryogenic 
equipment,  developed  a  series  of  dielectric  coatings  called  Laketites  for 
superconducting  wires  that  have  unusually  large  specific  heat  and/or  thermal 
conductivity  at  the  low  temperatures  that  are  characteristic  of  superconduc¬ 
tivity.  Simple  calculations  indicated  that  these  materials  would  assist  in 
solving  the  heating  problem  by  either  rapidly  carrying  away  the  heat  with 
high  thermal  conductivity  or  by  absorbing  heat  by  use  of  high  specific  heat 
or  by  some  combination  of  these  properties.  A  method  was  needed  to  choose  the 
best  material  in  this  series  for  further  testing  and  development  and  to  choose 
the  correct  thickness  of  the  coating  for  best  results.  Experimental  testing 
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of  different  coatings  and  thicknesses  would  require  testing  numerous  samples. 
Work  in  the  summer  of  19802  revealed  that  there  was  no  simple  experimental 
test  with  readily  available  equipment  that  could  screen  a  large  number  of 
samples.  As  an  alternative,  this  project  was  established  to  develop  a  com¬ 
puter  simulation  of  transient  heat  transfer  using  the  Air  Force-sponsored 
SCEPTRE  computer  code.  This  computer  program  was  then  used  to  evaluate 
different  coating  materials  and  various  thicknesses  of  coatings.  Preliminary 
reports  on  the  early  stages  of  the  development  of  this  program  have  been 
given  at  scientific  meetings. The  final  results  of  this  project  will  be 
presented  to  the  Applied  Superconductivity  Conference  in  Knoxville,  TN  November 
30,  1982  and  will  be  published  in  the  IEEE  Journal,  MAGNETICS  in  the  spring 
of  1982. 

This  report  describes  the  final  version  of  the  computer  program,  presents 
the  results  of  the  analysis  for  the  three  major  types  of  Laketite  coatings, 
and  gives  conclusions  about  their  feasibility  and  recommendations  for  future 
work. 


i,^r- 
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1 1 1  Theory 

The  basic  approach  of  this  project  is  to  use  the  powerful  and  versatile 
SCEPTRE  computer  code  to  simulate  transient  heat  transfer  processes.  The 
SCEPTRE  code  was  written  to  solve  fast  transient  electrical  problems.  To  solve 
heat  transfer  problems  requires  translating  the  thermal  process  into  electrical 
terms. 

Heat  transfer  in  a  solid  is  readily  modeled^  electrically,  as  a  series  of 
resistors  and  capacitors  with  the  resistors  representing  the  thermal  resistance 
due  to  a  finite  thermal  conductivity  and  the  capacitor  representing  the  heat 
capacity  of  the  object.  This  analogy  is  mathematically  exact,  i.e.  the  equa¬ 
tions  have  exactly  the  same  mathematical  form.  This  exactness  represents  the 
basic  similarity  between  the  two  processes  both  of  which  involve  transport  by 
particles--electrons  for  electricity,  phonons  for  heat.  Two  possible  problems 
can  arise  from  this  approach.  1)  A  basic  assumption  of  this  analysis  is  that 
the  mate-ials  are  isotropic--the  same  in  all  directions.  Some  materials 
notably  composites  are  not  isotropic.  In  the  absence  of  contrary  evidence,6 
it  is  assumed  that  the  Laketites  and  all  other  materials  involved  in  this 
project  are  isotropic.  2)  In  the  electrical  case  the  flow  of  current  is 
nearly  always  longitudinal  (the  long  direction)  which  usually  permits  the 
assumption  that  the  cross-sectional  area  perpendicular  to  the  flow  of  the 
current  is  constant.  The  flow  of  heat  is  often  radial  (the  short  direction): 

If  the  material  is  curved,  e.g.  a  wire,  the  cross  sectional  area  is  no  longer 
constant.  In  this  project,  the  simplest  possible  geometry--a  square  thin 
flat  plate  was  chosen.  Such  a  configuration  is  the  easiest  to  produce  ex¬ 
perimentally.  Thus  the  largest  amount  of  experimental  heat  transfer  data 
comes  from  thin  flat  plates.  Selecting  a  flat  plate  as  the  basic  configura¬ 
tion  allows  comparison  with  this  existing  data,  and  ease  in  setting  up  con¬ 
firming  experiments.  However,  most  actual  superconducting  projects  use 
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curved  wires.  The  program  can  be  readily  modified  to  solve  the  problem 
of  a  curved  wire.  Rather  than  use  actual  superconducting  material  as  the 
material  that  is  coated,  this  simulation  assumes  that  the  coating  is  applied 
to  copper.  All  superconductors  of  interest  are  imbedded  in  a  copper  or 
bronze  (copper  alloy)  matrix.  Since  copper  is  more  throughly  thermally 
characterized  than  are  superconductors  or  bronze,  copper  is  used  to  simulate 
the  superconductor. 

At  the  interface  between  the  material  and  the  liquid  helium,  there  is 
an  additional  boundary  resistance  commonly  called  the  kapitza  resistance. 

At  first  it  was  thought  that  this  resistance  only  appeared  between  materials 
and  superfluid  helium.  However,  experimental^  evidence  has  shown  that  this 
boundary  thermal  resistance  also  is  present  for  non-superfluid  helium.  Data 
between  materials  and  normal,  i.e.  nonsuperfluid,  helium  is  extremely  limited, 
but  indicates  that  the  resistance  is  of  the  same  order  of  magnitude  for  both 
normal  and  superfluid  helium.  To  account  for  the  kapitza  resistance  at  the 
liquid  hel ium--sol id  interface,  the  algebraic  form 

Rk  =  C/(T3(1+1.5AT/T+(AT/T)2  +  .25(AT/T)3))  (1) 

was  used.  This  is  derived  from  the  basic  theory  of  this  effect  while  allowing 

O 

the  assumption  that  AT  could  become  large  relative  to  T.  There  is  a  sound 
physical  theory  behind  this  form.  The  basic  theory  assumes  that  the  thermal 
resistance  at  the  liquid-solid  interface  arises  from  scattering  due  to  the 
difference  in  phonon  frequencies  in  the  material  and  the  helium,  i.e.  the 
acoustic  mismatch.  This  theory  predicts  that  the  principle  temperature 
dependence  of  the  resistance  should  be  1/T3  which  is  approximately  experi¬ 
mentally  correct.  It  furthermore  predicts  that  temperature  dependence  is 
independent  of  the  surface  condition  of  the  material  but  that  the  actual 

Q 

number  depends  on  surface  conditiors  which  is  in  agreement  with  observations. 

While  the  theory  must  b(  -~'Mflp  vo  achieve  detailed  agreement  with  experiment. 
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it  is  obviously  basically  correct.  Since  nearly  the  same  conditions  of 
acoustic  mismatch  hold  in  normal  helium  as  in  superfluid,  the  same  basic 
formula  should  be  correct.  In  superfluid  helium  the  kapitza  resistance 
dominates  heat  transfer  so  it  is  easily  measurable.  For  normal  helium  the 
kapitza  resistance  is  small  compared  with  other  effects,  thus  it  is  difficult 
to  measure.  Assuming  the  same  form  of  the  equation  and  coefficients  for  a 
copper-normal  helium  boundary  as  for  copper-superfluid  helium  boundary  seems 
to  be  the  most  reasonable  procedure.  The  acoustic  properties  of  the  Laketite 
have  never  been  measured,  as  they  are  new  material s.  Other  dielectries  show 
a  T“2  dependence  for  kapitza  resistance.^  However,  the  other  thermal  pro¬ 
perties  of  the  Laketites  are  similar  to  lead,  a  metallic  substance  with  a 
kapitza  resistance  similar  to  copper,  so  it  was  assumed  that  the  Laketites 
have  the  same  boundary  resistance  as  copper.  Nevertheless,  there  is  the 
possibility  that  this  assumption  may  not  be  correct  and  may  effect  the  results 
of  this  program.  If  the  actual  boundary  resistance  were  higher  for  the 
Laketites  than  the  program  predicts,  the  inner  surface  of  the  coating  under 
consideration  would  be  hotter  than  predicted  if  everything  else  remains  the 
same.  This  would  tend  to  make  both  the  superconductor  and  the  coating  become 
hotter  for  the  same  amount  of  heat  energy  input.  Conversely  if  the  actual 
kapitza  resistance  were  lower  than  the  predicted,  the  superconductor  would 
be  cooler  than  predicted.  However,  there  is  some  evidence11  that  hotter 
surface  temperature  due  to  thin  insulating  coating  e.g.  ice,  result  in  im¬ 
proved  heat  transfer;  because  the  insulator  promotes  the  local  superheating 
necessary  for  nucleate  boiling.  Simulating  this  problem  by  varying  the 
constant  in  equation  7  by  a  factor  of  six  changes  some  numbers,  but  does 
not  alter  the  general  form  of  the  results  or  any  of  the  basic  conclusions. 

Experimental  work  has  shown  that  there  is  a  kapitza  resistance  even 

12 

between  two  solids  such  as  epoxy  and  metal.  That  effect  was  visible  only 
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at  very  low  temperatures--^. IK,  well  below  the  temperature  of  boiling  helium-- 
M.2K.  This  solid  to  solid  contact  resistance  should  be  negligible  in  the 
higher  temperature  range,  4  to  10K,  of  this  report  unless  the  acoustic  pro¬ 
perties  of  Laketites  differ  substantially  from  that  of  copper.  Thus  any  such 
effect  was  omitted. 

One  should  note  that  the  technique  for  applying  the  SCI  compounds  that 
seems  most  promising  uses  a  composite  glass-Laketite  coating.  In  that  mixture 
the  thermal  conductivity  appears  to  be  dominated  by  the  kapitza  resistance 
between  the  glass  and  the  Laketite  coating;^  indicating  the  possible  complexity 
of  the  heat  transfer  situation.  In  the  absence  of  experimental  evidence, 
the  simplest  assumption-minimal  boundary  resistance--0  in  the  program—between 
the  superconductor  and  the  coating  seems  to  be  the  most  reasonable  approach. 
However  the  program  is  so  structured  that  any  copper-Laketi te  thermal  resis¬ 
tance  could  be  inserted  in  the  computer  simulation  by  adding  a  resistor 
between  the  copper  and  the  Laketites. 

The  addition  of  heat  to  the  wire-coating  system  is  simulated  by  an  ideal- 
no  resistance-current  source  distributed  throughout  the  length  of  the  sample. 

As  explained  in  the  introduction,  pulsed  operation  tends  to  produce  eddy 
current  heating  throughout  the  length  of  the  superconducting  wire.  Dielectric 
coatings  would  not  have  eddy  current  heating.  Neither  is  the  helium  subject 
to  heating  by  changing  magnetic  fields,  so  all  of  the  heat  source  is  in  the 
superconductor.  Eddy  current  heating  for  relatively  thin  matrixes  follows 
the  changing  magnetic  fields  with  a  time  constant  of  less  than  a  microsecond. 
Thus  for  times  longer  than  a  few  microseconds--a  square  wave  pulse  of  duration 
equal  to  the  total  switching  time  of  the  magnet  is  a  good  approximation--so 
it  was  used.  A  square  wave  pulse  assumes  linear  variation  of  the  magnetic 
field  during  switching  which  is  the  usual  approximation.  However  the  program 
is  constructed  so  that  any  form  of  heating  pulse  can  be  used. 

To  complete  the  equivalent  circuit  for  the  solid,  three  small  capacitors 
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were  inserted,  one  on  the  inside  edge  of  the  outermost  surface  of  the  copper, 
one  on  the  inside  edge  of  the  outermost  surface  of  the  coating  and  one  at  the 
innermost  edge  of  the  liquid  helium.  These  served  as  near  zero  mass  thermometers, 
at  those  three  points  since  the  program  would  automatically  compute  the  voltage 
across  these  capicitors--equivalent  to  temperature--at  all  times.  As  a  test 
of  their  effect  the  value  of  these  capacitors  was  increased  by  a  factor  of 
ten,  without  affecting  the  numerical  results,  at  the  one  tenth  of  one  per 
cent  level,  thus  allowing  the  conclusion  that  the  capacitances  of  these 
capacitors  had  no  effect  on  the  results. 

The  heat  transfer  to  the  liquid  helium  was  simulated  as  two  ideal  nega¬ 
tive  (removal  of  heat)  current  sources.  Only  a  thin  boundary  layer  of  helium 
would  be  heated,  but  the  amount  of  heat  transferred  into  and  out  of  the 
boundary  layer  was  varied  with  different  conditions.  The  value  of  the  heat 
transfer  from  the  solid  into  the  boundary  layer  was  controlled  by  a  lengthy 
subprogram  which  calculated  the  appropriate  amount  of  heat  transfer  based  on 
the  temperature  of  the  helium,  the  amount  of  energy  in  the  boundary  layer  of 
the  helium,  the  rate  of  flow  of  energy  into  the  helium,  the  rate  of  flow  of 
heat  from  the  boundary  layer  to  the  main  helium  bath,  the  elapsed  time,  and 
the  history  of  the  heat  transfer. 

Since  this  is  a  transient  heat  transfer  problem,  it  is  nearly  impossible 
to  add  enough  energy  to  heat  the  main  helium  bath  significantly.  Even  at  a 
phenomenal  50W/cm2  rate  for  a  100  cm2  surface,  a  .001  second  pulse  only  adds 
.05  J/cm2  to  the  bath,  sufficient  to  raise  one  liter  of  helium  only  .002K. 

The  main  bath  was  simulated  by  a  constant  ideal  voltage  source,  i.e. 
it  could  absorb  any  amount  of  heat  without  changing  temperature.  This 
temperature  was  fixed  at  4. OK  rather  than  the  normal  boiling  temperature  of 
4.21K  to  approximate  operation  under  slightly  reduced  pressure  as  would  be 
likely  in  an  aircraft  or  spaceship.  The  topological  properties  of  the 
SCEPTRE  code  requires  a  resistor  to  shunt  a  voltage  source.  The  value  of 
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this  "dummy"  resistor  was  varied  by  two  orders  of  magnitude  in  both  direc¬ 
tions  without  detectably  affecting  the  results  of  the  computer  program,  so 
its  presence  seems  unimportant.  The  complete  circuit  diagram  is  shown  in 
Figure  1. 

The  most  complex  problem  in  this  project  is  determining  the  amount  of 

heat  transfer  from  the  solid  into  the  boundary  layer  of  the  helium.  In 

steady-state  processes  there  are  four  distinct  modes  of  heat  transfer-- 

conduction  and/or  convection,  nucleate  boiling,  transition  boiling,  and 
5  9 

film  boiling.  ’  Transient  conditions  further  modify  some  of  these  modes. 

The  conduction-convection  region  is  the  most  heavily  modified  by  tran¬ 
sient  conditions.  In  the  steady  state  low  heat  flow  is  dominated  by  convec¬ 
tion.  Convective  heat  transfer  is  difficult  to  simulate  since  it  strongly 
depends  on  the  geometry  of  the  situation;  e.g.  the  orientation  of  the  heater, 
the  size  of  the  channels,  the  depth  of  the  fluid,  etc.  However,  for  short 
time  periods  conductive  effects  dominate.  Convection  is  driven  by  density 
differences  in  a  gravitational  field.  A  maximum  for  the  acceleration  of  a 
parcel  of  fluid  next  to  the  solid  is  of  the  order  of  .Olg  where  g  is  the 
acceleration  of  gravity.  Assuming  continuing  constant  acceleration  of  .Olg, 
the  maximum  distance  that  hot  fluid  could  travel  in  .001  seconds  would  be 
4.9  x  10'6cm.  This  maximum  distance  of  travel  also  represents  the  greatest 
reasonable  depth  of  fluid  that  can  be  influenced  by  convective  heat  transfer. 

The  maximum  distance  that  is  influenced  by  conductive  heat  transfer  is  approx¬ 
imately  given  by  at  where  a  is’the  thermal  diffusivity  K/pcp.  For  helium  at 
4K,  a=3.2  x  lO^cmVsec.  So  for  a  time  of  .001  seconds  the  conductive  heat 
penetration  depth  is  about  5.66  x  10_4cm,  much  farther  than  the  convective 
distance.  At  about  .05  second  both  distances  are  of  the  same  order  of  magnitude. 
On  the  other  hand  for  very  short  times  the  distance  that  heat  penetrates  begins 
to  approach  the  depth  of  the  surface  roughness  which  is  usually  of  the  order 
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FIGURE  1 


Electrical  circuit  no del  of  liquid  helium 
in  contact  with  a  heated  surface. 
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Figure  lb.  Detail  model  of  Coating  and  Superconductor. 
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of  10“5cm  which  invalidates  an  analysis  based  on  flat  plates  and  flat 
boundary  layers.  Thus  between  some  very  short  time  and  about  10-2  seconds, 
conductive  heat  transfer  is  the  dominant  mode  and  convection  can  be  ignored. 

Before  considering  the  exact  mechanism  for  conductive  heat  transfer,  the 
thickness  of  the  boundary  layer  must  be  determined.  The  lower  limit  of  this 
size  can  be  set  two  ways.  Even  the  highly  polished  surfaces  have  irregular¬ 
ities  of  the  order  of  10"5cm,  i.e.  10000A.  Unless  the  boundary  layer  is 
substantially  thicker  than  the  irregularities,  a  model  of  a  flat  boundary 
layer  on  a  flat  surface  is  meaningless.  Also,  it  is  desirable  to  have  the 
same  boundary  layer  serve  for  nucleate  boiling  as  well  as  conduction.  The 
minimum  diameter  for  bubbles  in  boiling  in  Helium  is  10“ 5  to  10“Gcm  depending 
on  the  theory  used.  The  maximum  size  of  the  boundary  layer  is  set  by  thermal 
diffusion  constraints.  Because  a  is  small  for  helium,  the  maximum  thermal 

diffusion  ranges  from  about  5  x  10_4cm  for  one  millisecond  to  about  1  x  10_3cm 
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for  10  milliseconds.  The  Japanese  work  showing  that  grooves  of  less  than 
1mm  width  aid  heat  transfer  can  easily  be  understood  as  an  increase  in  effec¬ 
tive  area  if  the  boundary  layer  is  substantially  less  than  one  half  the  groove 
width.  If  the  boundary  layer  is  thicker  than  one  half  the  groove  width,  the 
boundary  layers  from  opposite  sides  of  the  groove  touch  and  there  is  no  in¬ 
crease  in  effective  area.  Even  a  close  approach  of  the  boundary  layers  would 
inhibit  heat  transfer  in  the  same  manner  that  narrow  channels  inhibit  boiling 
heat  transfer.  The  evidence  argues  for  a  maximum  depth  of  the  boundary  layer 
of  10~3cm.  Splitting  the  difference  yields  a  depth  of  2  x  10'4cm  for  the 
boundary  layer.  As  a  test  of  the  sensitivity  of  the  program  to  this  parameter, 
the  depth  of  the  boundary  layer  was  both  increased^and  decreased  from  this 
value  by  a  factor  of  5  without  affecting  the  results. 

Having  determined  the  depth  of  the  boundary  layer,  the  procedure  was  to 
construct  a  temperature  profile  for  the  boundary  layer.  Heat  transfer  out 
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of  the  boundary  layer  was  given  by 
Q/A  -  K  dT/dx 


(2) 


with  the  derivative  being  evaluated  at  90%  of  the  depth  of  the  boundary 
layer  to  avoid  mathematical  singularities  at  the  exact  boundary.  A  temperature 
profile  of  the  form 


AT’(x)  =  AT..  ,(l-x/D)y 


total 

where  AT'(x)  is  the  rise  in  temperature  above  the  original  bath  at  a  dis¬ 
tance  x  away  from  the  solid  surface,  AT^^-j  is  the  total  difference  in 
temperature  between  the  solid  and  the  bath,  D  is  the  thickness  of  the  boundary 
layer  and  y  is  an  exponent  to  be  determined,  y  is  determined  by  requiring 


(3) 


that  the 

0D  pCpAT'(x)  dx  =  E  (4) 

where  E  is  the  total  thermal  energy  present  in  the  boundary  layer,  p  is 
the  density  and  Cp  the  specific  heat  of  the  helium  in  the  boundary  layer. 

Since  E  is  known  (the  computer  is  very  efficient  at  keeping  track  of  all 
heat  inputs  and  outflows  from  the  boundary  layer  and  maintaining  a  running 
total  of  the  energy  contained  in  the  boundary  layer),  substituting  Equation 
3  in  Equation  4  determines  the  exponent  y.  With  y  determined,  Equation  3 
can  be  used  to  determine  T(x).  Differentiating  T(x)  and  using  Equation  2 
determines  the  heat  flux  out  of  the  boundary  layer.  Clearly  heat  flux  out 
of  the  boundary  layer  is  conductive  for  short  periods  of  time  since  convec¬ 
tive  motion  can't  reach  the  outer  edge  of  the  boundary  layer.  In  steady 
state  both  convection  and  conduction  form  linear  temperature  profiles  which 
is  the  limiting  case  for  this  procedure.  The  details  of  heat  transfer  inside 
the  boundary  layer  are  not  known.  The  use  of  energy  consideration  to  deter¬ 
mine  heat  flow  in  the  absence  of  detailed  knowledge  is  customary  procedure. 

A  polynomial  curve  is  the  simplest  equation  that  reduces  to  the  proper  steady 
state  limit.  As  Figure  2  shows,  the  temperature  profiles  generated  by  this 


Temperature  Profile  in  the  Thermal  Boundary 
Layer  for  Conduction  Heat  Transfer. 
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method  yield  maximal  heat  flow  into  the  boundary  layer  and  very  little  out¬ 
flow  of  heat  for  short  periods  of  time  and  approach  steady  state  for  longer 
periods  of  time.  Thus  this  method  of  calculating  heat  flux  out  of  the  boundary 
layer  seems  physically  reasonable.  Unfortunately  this  method  can  not  handle  the 
situation  when  the  surface  is  cooling,  i.e.  recovery. 

As  the  temperature  of  the  solid  surface  increases  and  the  energy  con¬ 
tained  in  the  boundary  layer  increases,  eventually  nucleate  boiling  commences. 

The  traditional  view  of  the  onset  of  nucleate  boiling  is  that  gas  trapped  in 
minute  cavities  on  the  surface  expands  to  form  bubbles,^  a  process  called 
heterogeneous  nucleation.  The  onset  of  nucleation  requires  that  temperature 
T^  for  the  liquid  adjacent  to  the  wall  be  greater  than  the  boiling  temperature 
of  the  fluid  T$at.  Most  sources  use  T  but  that  ignores  the  kapitza 
resistance.  In  addition,  if  the  bubble  is  to  expand  and  grow  there  must  be 
sufficient  energy  present  for  some  of  the  liquid  around  the  bubble  to  vaporize. 
Since  a  single  bubble  does  not  affect  the  heat  transfer  properties  by  a 
measurable  amount,  exact  criteria  for  the  onset  of  nucleate  boiling  are  hard 
to  develop.  Such  criteria  are  necessary  for  a  computer,  which  has  trouble 
handling  a  situation  where  the  mode  of  computation  is  an  indefinite  mixture 
of  two  situations. 

For  hetrogenous  nucleation  a  criteria  for  the  beginning  of  nucleate 
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boiling  can  be  derived  from  the  equation  by  Han  &  Griffith  for  the  critical 
radius  of  nucleation  sites. 

rc  "  <TL  '  Tsat)  !+  1-  12  <TL  '  Tbath>  Tsat  a  (5) 

-  Tbath>  <tl  •  Tsat>J 

The  critical  radii  first  became  real  when  the  term  under  the  square  root 
equals  or  exceeds  zero.  For  the  case  where  Tbatb  =  Tsat  the  beginning  of 
nucleate  boiling  becomes 
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<TL  -  Tbath>  =  — Wl 


(6) 


pv  hfg 

This  can  be  related  to  energy  by  multiplying  by  Cp,  giving 


-NB 


12pL  Cp  Tbath  0 
pv  hfg 


(7) 


For  saturated  liquid  helium  at  4K  the  energy  for  the  beginning  of  nucleate 
boiling  is  E^g  =  1.04  x  10~2  Jm-2.  To  allow  for  the  small  range  of  radii 
necessary  to  actually  produce  boiling,  a  value  of  2  x  10"2Jm“2  was  used  for 


^BNB  in  this  Pr°9ram-  This  value  gave  results  that  agreed  very  well  with 
experimental  data.  The  computer  added  up  the  total  energy  in  the  boundary 
layer;  when  this  figure  exceeded  Eg^g  nucleate  boiling  was  assumed  to  occur. 
However,  experimentally  nucleate  boiling  does  not  begin  with  a  sharp  transi¬ 
tion,  but  rather  with  a  mixture  of  both  nucleate  boiling  and  conductive 
heat  transfer.  The  smooth  joining  of  the  conductive  and  the  nucleate  regime 
was  achieved  as  follows.  For  energies  in  the  boundary  layer  of  less  than  four 
times  the  minimum  energy  required  for  the  initiation  of  nucleate  boiling,  the 
heat  transfer  is  calculated  as  a  linear  combination  of  conductive  and  nucleate 
boiling  heat  transfer. 

Very  recent  experimental  work*6’*7  has  suggested  that  hetrogeneous 
nucleation  is  incorrect  for  helium  because  of  the  absence  of  trapped  gases. 

If  hetrogeneous  nucleation  does  not  occur,  nucleate  boiling  still  exists 
commencing  with  a  process  called  homogenous  nucleation  which  occurs  when  local 
temperature  and  energy  flucuations  create  bubbles  in  the  liquid.  The  onset 
of  nucleate  boiling  is  very  sharp  and  occurs  solely  as  a  function  of  tempera¬ 
ture  for  a  liquid  under  fixed  pressure.  The  difference  in  the  conditions  for 
onset  of  nucleate  boiling  between  the  two  theories  is  small  (barely  detectable 
experimentally).  Nevertheless  future  development  of  this  computer  program  should 
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include  a  changing  of  the  criteria  for  the  onset  of  nucleate  boiling  to  the 
homogenous  nucleation  case  since  both  the  physical  arguments  and  the  experi¬ 
mental  evidence  for  the  existence  homogenous  nucleation  in  helium  appear 
compelling. 

Once  nucleate  boiling  becomes  established  it  dominates  heat  transfer 
because  the  relation  between  q  and  AT  becomes  much  greater  than  the  linear 
relationship  for  conductive  heat  transfer.  Nucleate  boiling  once  well  estab¬ 
lished  is  relatively  independent  of  the  material  or  the  condition  of  the 
surface  of  the  solid.  Correlations  such  as  the  one  Kutateladze  developed 


based  solely  on  the  hydrodynamic  properties  of  the  fluid  agree  fairly  well 

with  a  wide  range  of  experimental  data.  Using  the  available  experimental 
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data  and  guided  by  the  general  temperature  variation  of  the  quantities  in 
the  hydrodynamic  correlations,  a  form 

»  ■  B<TW-Tsat>!  (' 

was  used  as  the  basic  formula  for  heat  transfer  in  the  nucleate  boiling  region. 

The  constant  is  evaluated  by  requiring  that  a  q  vs  T  plot  gives  the 
correct  results  at  the  peak  of  the  nucleate  boiling  curve.  Because  of  the 
dramatic  changes  in  heat  transfer  at  the  peak  of  nucleate  boiling,  this  point 
is  easily  determined  experimentally.  Data  for  various  copper  systems  similar 
in  configuration  to  the  arrangement  in  this  program  cluster  around  the  values 
used  in  this  program  which  are  q/A  =  .6  W/cm2  and  -  T$at  =  .5K.  The 
higher  values  reported  in  some  work  are  somewhat  suspect  since  no  allowance 
was  made  for  kapitza  resistance. 

At  the  peak  heat  flux  for  nucleate  boiling,  the  bubbles  merge  in  an 
unstable  fashion  to  form  a  continuous  layer  of  gas.  Because  the  thermal 
conductivity  of  the  vapor  is  much  lower  than  that  of  the  liquid--a  factor 
of  five  for  helium--a  large  resistance  to  heat  flux  builds  up.  With  the 
build-up  of  this  resistance  the  amount  of  heat  transported  across  the 
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boundary  layer  into  the  helium  bath  drops  rapidly.  If  the  heat  flux  remains 
constant,  the  temperature  of  the  sample  rises  rapidly,  a  phenomenum  called 
burnout.  Due  to  the  rapidity  of  the  temperature  rise,  typically  a  few  micro¬ 
seconds,  the  exact  form  of  the  heat  transfer  in  this  region  is  unknown.  From 
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measurements  with  the  surface  temperature  carefully  controlled  it  is  known 
that  the  heat  flux  for  copper  would  drop  to  about  q/A  =  .2W/cm2  at  a  temperature 
difference  -  T  =  6.8K.  Until  recently  little  study  of  this  region  has 
been  made,  so  that  a  linear  fit  between  the  peak  flux  and  the  minimum  flux 
seems  to  be  the  best  approach  for  the  q  vs  AT  curve.  This  region,  which  is 
known  as  transition  boiling,  was  modeled  as  a  combination  of  nucleate  and  film 
boiling  heat  transfer.  It  was  assumed  that  transition  boiling  occurred  be¬ 
tween  the  peak  nucleate  boiling  flux  and  the  establishment  of  film  boiling. 
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Following  the  theory  of  Kalinin  the  fraction  of  nucleate  boiling  contributing 
to  the  heat  transfer  is  calculated  from 


f  =  (1  -AT*)  (9) 

where 

AT*  =  ((FBL(VTsat)3)‘25  '  (CT1}-25  '  ((ct2)*25  "  (Cn)-25)  (10) 
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where  CT1  and  CT2  are  constants.  The  data  of  Steward  for  carbon  surfaces 

were  used  to  determine  C-^  =  1.  and  =  15  x  10‘‘  J  K3/m2  respectively  for 

the  incipience  and  the  completion  of  transition  boiling. 

To  avoid  the  problems  of  the  transition  region,  several  attempts  have 

been  made  to  define  a  criterion  for  transition  from  nucleate  boiling  to  film 
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boiling  during  transient  heat  transfer.  None  of  the  approaches  fit  the 

data  very  well. 

The  most  promising  approach  seems  to  focus  on  the  energy  contained  in  the 
boundary  layer.  For  short,  transient  heat  pulses,  less  than  .01  seconds  in 
duration,  the  bubbles  from  nucleate  boiling  do  not  have  sufficient  time  to  escape. 
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Since  the  energy  cannot  escape  from  the  boundary  layer,  a  simple  criterion 

for  the  start  of  film  boiling  is  that  the  energy  in  the  boundary  layer  be 

equal  to  or  exceed  the  amount  of  energy  required  to  vaporize  all  of  the 

liquid  in  the  boundary  layer.  It  produced  a  formula  E  Tm  =  a  constant  with 

m  between  2  and  2.5  depending  on  the  exact  assumption  of  the  behavior  of  the 

boundary  layer.  This  approach  has  limited  agreement  with  the  available 
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experimental  data  but  not  sufficient  agreement  to  use  in  the  computer 
program.  Instead  an  empirical  formula 

E8L  (TL  '  Tbath)!  =  constant  (11) 

with  the  constant  being  5  x  104  Jm-2  was  used,  which  agrees  somewhat  better 
with  the  experimental  data. 

Steady  state  film  boiling  heat  transfer  is  best  modeled  by  the  cor- 
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relation  of  Breen  &  Westwater  .  This  must  be  modified  to  account  for  tran- 
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sient  effects  as  pointed  out  by  Iwasa.  To  model  these  transient  conditions 
the  following  form  was  used 

qFB  -  (200(Tl  -  Tbath)>.>)  (  (Efbs$  ♦  Ebl)/  2Ebl)  (12) 

since  energy  variations  are  more  easily  calculated  in  the  SCEPTRE  format 
than  are  temperature  variations. 

Actual  techniques  of  the  use  of  the  SCEPTRE  code  are  contained  in  the 
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instruction  manual  and  the  book  on  the  code.  An  example  of  an  actual 
SCEPTRE  program  is  given  in  Appendix  I  and  the  critical  Fortran  subroutine 
that  actually  calculates  the  value  of  the  boiling  heat  transfer  is  given  in 
Appendix  II. 
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IV  Results 

With  the  computer  program  operational,  the  procedure  was  to  run  pro¬ 
grams  using  different  thicknesses  of  coatings,  different  compositions  of 
coatings,  different  levels  of  heating  and  different  times.  Based  on  the 
manufacturer's  report  on  the  Laketites,^  there  are  three  broad  classes 
of  these  materials,  SC-1,  SC-2,  SC-3  with  several  different  variants  in 
each  class.  The  most  promising  compositions  are  SC-1C,  SC-2B,  and  SC-3A 
so  these  are  the  ones  that  were  studied.  For  comparison  purposes  the  same 
computer  runs  were  made  on  bare  copper  (no  coating)  and  on  copper  that  was 
coated  with  GE  7031  varnish.  GE  7031  is  typical  of  the  usual  dielectric 
materials  used  to  insulate  wires  at  low  temperatures.  Because  of  its  ease 

of  application,  it  is  the  most  popular  such  material  and  therefore  has  been 
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the  subject  of  many  investigations  to  accurately  determine  its  low  temp¬ 
erature  properties.  In  the  case  of  the  SC-1  materials,  GE  7031  is  mixed 
with  the  coating  material  in  the  application  process.  Thus  the  logical  choice 
as  the  comparison  material  representing  typical  existing  dielectric  coatings 
at  low  temperatures  is  GE  7031.  The  superconductor  is  represented  by  a  slab 
of  copper  .005"  (5mils)  =  .0125cm  thick.  This  is  a  typical  thickness  for 
the  copper  or  bronze  matrix  around  the  superconducting  filaments  in  the  normal 
construction  of  filamentary  superconducting  cable.  Coating  thicknesses  of 
1,  2,  4,  10,  and  20  mils  were  used  (.0025,  .005,  .01,  .025,  .05cm).  Mechanical 
constraints  in  the  process  of  applying  coatings  to  wires  limit  the  coating 
thickness  to  this  range  of  values.  Superconductors  might  face  a  wide  variety 
of  heating  levels  from  gentle  heating  for  slow  switching  and  low  magnetic 
fields  to  severe  heating  for  rapid  switching  and  high  magnetic  fields.  This 
was  simulated  by  varying  the  input  power  level  of  heating  in  the  copper.  The 
heating  levels  chosen  were  0.05,  0.1,  0.2,  0.5,  1.0,  2.0,  3.0,  4.0,  5.0,  10.0, 
20.0  and  50.0  W/cm2.  The  clustering  of  values  in  the  1-5  W/cm2  range  reflected 
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crude  calculations  that  this  was  the  actual  heating  range  for  some  specific 
Air  Force  projects  that  are  in  the  design  stage.  All  values  are  in  W/cm2, 
the  usual  scientific  convention  for  removing  the  effect  of  varying  areas  in 
different  experimental  arrangements.  All  heating  was  assumed  to  be  a  constant 
heating  rate,  with  zero  turn-on  time,  i.e.  a  square  wave  pulse  with  a  zero 
turn-off  time  at  the  end  of  the  simulation.  This  is  a  good  approximation 
or  real  eddy  current  heating,  but  real  experiments  have  a  finite  turn  on 
time,  usually  of  the  order  of  2  to  20  microseconds,  thus  a  comparison  of  these 
results  with  experiments  where  the  heating  is  less  than  10-'4  seconds  in  duration 
is  questionable  unless  a  correction  is  introduced  to  account  for  the  finite 
turn-on  time.  The  length  of  heating  and  the  total  time  for  the  computer 
simulation  of  the  behavior  were  identical  since  this  program  can  not  success¬ 
fully  model  the  problem  of  recovery--the  cooling  down  after  the  end  of  a 
heating  pulse.  Three  different  times  were  used--0.001  seconds,  0.01  seconds 
and  0.1  seconds.  For  time  periods  substantially  shorter  than  1  millisecond 
the  actual  surface  roughness  (usually  an  unknown  variable)  must  be  accounted 
for  as  the  approximation  of  a  flat  surface  and  a  flat  boundary  layer  breaks 
down.  For  times  longer  than  100  milliseconds,  the  basic  assumption  of  a 
constant  thickness  boundary  layer  breaks  down  because  of  convective  motion. 

Thus  these  three  times  span  the  time  range  that  is  normally  meant  by  the  term 
transient  heat  transfer. 

As  explained  in  the  theory  section,  it  was  assumed  that  heating  was 
uniform  throughout  the  copper  due  to  eddy  current  heating  from  changing  mag¬ 
netic  fields;  thus  all  the  heat  flow  was  assumed  to  be  radial.  The  important 
case  of  localized  heating,  'hot  spot1,  where  both  radial  and  longitudinal 
heat  flow  would  be  present  was  not  considered.  The  'hot  spot'  problem 
requires  a  two  dimensional  solution  to  the  heat  equations  which  is  beyond 
the  scope  of  this  project.  However,  another  project  which  would  investigate 


this  problem  would  have  considerable  merit,  since  the  results  of  this  pro¬ 
ject  indicated  that  these  materials  are  potentially  very  useful  for  distri¬ 
buted  heating  cases  and  would  appear  to  also  provide  considerable  protection 
against  localized  'hot  spots'.  If  the  insulation  coating  could  help  provide 
thermal  protection  against  'Hot  spots',  the  bulk  of  the  superconducting  device 
could  be  reduced  and  the  reliability  of  these  devices  enhanced. 

The  three  classes  of  Laketites  materials  are  primarily  distinguished 
thermally  by  their  different  specific  heats  and  thermal  conductivity  at  low 
temperatures.  The  SC-1  materials  have  very  high  specific  heat,  in  the  5  to 
10K  range,  and  thermal  conductivity  that  is  relatively  low,  only  slightly 
higher  than  7031  and  other  currently  existing  dielectric  materials.  The 
SC-2  materials  have  relatively  high  specific  heat,  though  not  as  large  as 
the  SC-1  materials,  and  very  high  thermal  conductivity,  up  to  one  third  of 
that  of  metallic  copper  (copper  is  the  universal  choice  for  the  highest  levels 
of  thermal  conductivity  at  low  temperatures) .  The  SC-3  materials  have  extremely 
high  thermal  conductivity,  approaching  that  of  metallic  copper,  but  have 
relatively  low  specific  heat. 

The  results  of  this  project  clearly  showed  that  for  the  distributed 
continuous  heating  case  thermal  conductivity  is  the  most  important  parameter. 
Figure  3  gives  an  example  with  a  superimposed  plot  of  the  change  of  tempera¬ 
ture  of  the  copper,  i.e.  the  superconducting  wire,  versus  time.  Each  coating 
is  two  mils  thick,  the  power  level  is  1  W/cm2  and  the  time  is  .01  seconds. 

The  wire  coated  with  SC-1C  reaches  a  final  rise  in  temperature  of  9.8K,  i.e. 
a  final  temperature  of  13. 8K,  which  exceeds  the  transition  temperature  of 
most  superconductors  including  all  of  the  commonly  used  superconductors  for 
magnets.  In  addition,  the  temperature  of  the  superconductor  rises  rapidly, 
reaching  a  value  of  5K  above  bath  temperature  in  1.3  milliseconds.  Five  K 


TIME  (ms) 


Pig. 3.  The  temperature  of  the  superconductor  (TWIRE) 
as  a  function  of  time  for  an  input  power  of 
lVI/cnr  for  wire  coated  with  different  thicknesses 
of  SC-1C  coating. 
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brings  the  temperature  to  9K,  which  is  the  absolute  0  field  limit  for 
Nb^Sn  wire.  In  a  real  magnet  such  a  wire  would  go  normal  due  to  the  residual 
magnetic  field.  Nb^Sn  is  the  highest  transition  temperature  superconducting 
magnet  wire,  so  a  temperature  rise  of  5K  is  certain  trouble.  In  contrast, 
wire  coated  with  SC-2B  has  warmed  only  .98K,  to  a  temperature  of  4.98K, 
after  the  full  10  milliseconds  of  heating--wel 1  within  the  limits  of  operation 
of  all  commonly  used  superconductors.  Running  the  program  for  a  longer  time 
reveals  that  SC-2B  caused  the  wire  to  heat  to  5K  only  after  22  milliseconds, 
over  an  order  of  magnitude  difference  from  SC-1C.  Clearly  there  is  a  tremen¬ 
dous  improvement  in  performance  for  SC-2B  vs  SC-1C.  Even  at  a  power  of 
0.5  W/crir  SC-1C  coated  wire  heats  beyond  the  crucial  5K  level  in  3.8  milli¬ 
seconds.  Even  at  power  levels  as  low  as  0.1  W/cm2  the  total  temperature  rise 
for  10  milliseconds  for  SC-1C  is  an  uncomfortably  high  2.7K. 

The  results  of  SC-1C  aid  7031  are  quite  similar  to  each  other,  as  are 
the  results  for  the  grouping  SC-2B,  SC-3A  and  bare  copper.  SC-1C  and  7031 
share  similar  low  thermal  conductivities  while  differing  sharply  in  specific 
heat.  SC-2B,  SC-3A  and  bare  copper  have  high  thermal  conductivity  out  to 
the  solid-liquid  interface  but  differ  substantially  in  specific  heat. 
Furthermore,  the  temperature  of  the  boundary  layer  of  helium  is  remarkably 
similar  (within  .05K)  for  all  coated  and  bare  surfaces  for  all  power  levels 
below  2W/cnr  and  within  10%  agreement  at  higher  power  levels.  The  tempera¬ 
ture  rise  of  wires  coated  by  SC-1C  and  GE  7031  seems  to  be  caused  by  large 
temperature  drops  required  to  transport  energy  across  low  thermal  conductivity 
materials. 

The  same  result,  that  thermal  conductivity  is  the  dominant  parameter, 
is  evident  by  comparing  results  with  the  same  material  but  differing  thick¬ 
nesses  of  coating.  The  total  heat  capacity  depends  on  the  product  of  the 
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volumetric  specific  heat  and  the  volume  of  the  material.  The  volume  of 
the  coating  is  directly  proportional  to  the  thickness  of  the  coating.  As 
the  thickness  of  the  coating  is  varied  from  1  to  20  mils,  the  volume  and  the 
total  heat  capacity  change  by  a  factor  of  twenty.  In  contrast  the  effect 
of  high  thermal  conductivity  is  lessened  by  additional  thickness  since  the 
temperature  change  across  a  material  is  directly  proportional  to  the  thick¬ 
ness  of  the  material.  Thus,  a  thicker  material  will  'soak  up'  more  heat 
as  it  becomes  thicker,  causing  the  superconductor  to  have  a  smaller  tempera¬ 
ture  rise;  but  the  increased  temperature  difference  required  for  thermal 
conduction  due  to  the  increased  thickness  causes  a  larger  temperature  rise 
of  the  superconductor.  If  specific  heat  is  the  dominant  factor,  the  super¬ 
conductor  will  become  much  cooler  for  a  given  level  as  the  thickness  is 
increased. 

As  is  shown  in  the  superimposed  plot  for  differing  thickness  of  SC-1C, 
this  high  specific  heat  material  does  have  a  drop  in  superconductor  tempera¬ 
ture  and  a  changing  rate  of  heating  with  increasing  thickness.  A  further 
example,  is  that  for  a  power  level  of  lW/cnr  for  a  1  ms  pulse,  the  maximum 
temperature  rise  of  the  superconductor  falls  from  5. OK  for  a  1  mil  coating 
of  SC-1C  to  2.55K  for  a  20  mil  coating.  However,  this  is  a  small  effect 
for  a  20  times  increase  in  volume,  particularly  since  under  the  same  condition 
SC-2B--a  medium  specific  heat  coating--causes  the  temperature  rise  of  the 
superconductor  to  decrease  from  0.74K  for  a  1  mil  coating  to  0.68K  for  a 
20  mil  coating  and  SC-3A--a  lower  specific  heat  material--shows  an  increase 
of  only  .OIK  in  the  temperature  rise  of  the  superconductor  from  0.74K  to  0.75K. 
If  specific  heat  were  the  key  variable,  a  20  mil  coating  of  an  SC-1  material 
should  give  improved  performance  over  a  1  mil  coating  of  SC-2  or  SC-3  material. 
This  does  not  happen.  The  high  thermal  conductivity  materials  outperform 
the  ones  with  low  thermal  conductivity  at  all  thicknesses.  Nevertheless, 


Fig. 4.  The  temperature  of  the  superconductor  as  a  ? 

function  of  time  for  an  input  power  of  117/ crrr 
for  wire  coated  with  different  thicknesses  of 
SC-1C  coating. 


27 


SC-1  materials  in  thick  coatings  are  definitely  superior  to  7031  whose 
performance  worsens  with  greater  thickness  due  to  its  very  low  specific 
heat  and  thermal  conductivity.  These  results  suggest  that  if  SC-1C  is 
the  material  of  choice  for  further  development  that  the  thickest  possible 
coatings  be  used  to  best  utilize  its  high  specific  heat. 

The  power  level  of  the  heating  might  affect  the  choice  of  coating  material. 
Table  I  presents  the  maximum  temperature  rise  for  the  superconductor  for  2  mil 
coating  for  1  ms  pulse  for  SC-1C,  SC-2B,  SC-3A  bare  copper,  and  7031. 

TABLE  I 

Maximum  temperature  rise  of  the  superconductor  versus  heater  power  for  a  1  ms 
pulse  and  a  coating  thickness  of  2  mil. 


Power 

W/cirr 

SC-1C 

K 

SC-2B 

K 

SC-3A 

K 

Bare  Cu 

K 

7031 

K 

.05 

0.8 

.30 

.32 

.32 

.68 

1 

1.45 

.35 

.37 

.36 

1.10 

.2 

2.25 

.42 

.42 

.42 

1.85 

.5 

3.45 

.52 

.52 

.52 

3.8 

1 

4.30 

.72 

.76 

.74 

6.0 

2 

7.1 

1.30 

1.35 

1.35 

8.8 

3 

9.2 

2.10 

2.6 

3.0 

11.0 

4 

U.O 

3.6 

6.2 

9.9 

12.5 

5 

12.0 

5.0 

8.0 

11.0 

14.0 

10 

17.0 

9.9 

11.5 

15.5 

18.0 

20 

22.5 

14.0 

15.1 

20.5 

24.0 

50 

31 

21.0 

22 

29 

32.0 

Clearly,  at  all  power  levels  the  high  thermal  conductivity  materials 
SC-2B,  SC-3A  outperform  the  low  thermal  conductivity  ones,  by  having  a  lower 
temperature  rise.  The  greater  specific  heat  of  SC-1C  begins  to  improve  its 
performance  relative  to  7031  at  higher  power  levels,  but  otherwise  specific 
heat  seems  to  have  little  effect  on  the  temperature  rise.  The  results  for  all 


29 


program  assumes  that  the  solid-helium  interface  is  being  substantially 
heated  and  has  difficulties  when  no  energy  is  reaching  the  interface.  Since 
a  is  a  key  in  the  treatment  of  isolated  hot  spots,  it  would  appear  that  the 
large  a  materials  SC-2B  and  SC-3A  would  also  give  increased  protection  against 
localized  heating. 

Time  could  also  be  a  factor.  Table  III  shows  the  maximum  temperature 
rise  for  2  mil  thick  coating  at  a  power  level  of  .5  W/cm2  for  different 
times  for  each  material. 


TABLE  III 


Time 

lxl0_3sec 

SC-1C 

K 

SC-2B 

K 

SC-3A 

K 

7031 

K 

1 

3.45 

.52 

.52 

3.8 

10 

6.5 

.52 

.52 

4.0 

100 

6.5 

.52 

N.A. 

4.1 

Note:  at  100ms  the  SCEPTRE  program  for  SC-3A  had  excessive  computation  time 
and  failed  to  reach  this  power  level  before  exceeding  the  pre-set 
computer  time  limits. 

The  relative  performance  of  the  coating  seems  to  be  unaffected  by  the 
length  of  the  heating  pulse.  Obviously  for  a  power  level  of  0.5W/cnr  the 
SC-2B  and  SC-3A  are  acceptable,  SC-1C  and  7031  are  not. 

The  computer  program  generated  far  more  data--over  10,000  pages  of 
printout  and  graphs--than  can  be  presented  in  this  report.  All  of  the  computer 
data,  both  numbers  and  graphs,  has  been  stored  at  Saginaw  Valley  State  College 
and  is  available  to  any  interested  party.  The  basic  SCEPTRE  programs  for 
each  coating  and  the  Fortran  subroutines  have  been  stored  on  magnetic  tape 
at  the  Michigan  State  University  Computer  center.  The  basic  SCEPTRE  language 
tapes  are  also  stored  there.  In  addition,  all  of  the  programming  information 
is  stored  at  Wright-Patterson  Air  Force  Base.  These  programs  can  be  reactivated 
on  very  short  notice  should  new  materials  be  developed  or  materials  other  than 
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other  coating  thicknesses  are  essentially  the  same.  Thicker  coatings  usually 
lead  to  lower  temperature  rises  but  the  high  thermal  conductivity  materials 
have  a  much  lower  maximum  temperature  rise  for  all  power  levels. 

Since  this  is  a  transient  process  one  might  except  that  the  correct  para¬ 
meter  for  dealing  with  it  is  the  thermal  diffusivity  a  =  K/pc ,  the  values  of 
which  are  given  in  Table  II  for  a  temperature  of  5K--roughly  the  mid-point 
of  the  heating  process. 

TABLE  II 

Values  of  a  for  different  materials  at  5K 


The  better  behavior  is  associated  with  the  materials  with  larger  thermal 
diffusivity,  as  expected.  But  SC-3A  and  SC-2B  have  nearly  identical  thermal 
response  despite  large  differences  in  thermal  diffusivity.  From  the  definition 
of  thermal  diffusivity  it  is  obvious  that  the  most  desirable  property  is  large 
thermal  conductivity  and  low  specific  heat.  Judging  from  the  computer  results 
this  is  indeed  the  correct  combination. 

Another  indication  of  the  importance  of  the  thermal  diffusivity  is  the 
unusual  results  for  very  thick  coatings  of  SC-1C.  The  surface  where  the  liquid 
helium  meets  the  coating  develops  small  temperature  oscellations,  including 
negative  temperature  rises--i.e.  temperature  falling  below  the  bath  tempera¬ 
ture--  for  thick  (10  &  20  mil)  SC-1C  coatins  at  low  power  levels.  The  oscil¬ 
lations  are  small  (.001K)  and  appear  to  be  caused  by  the  inability  of  the 
heat  pulse  to  reach  the  surface  due  to  the  low  value  of  a  for  SC-1C.  This 
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7031  need  to  be  tested  for  comparison  purposes.  Further  modifications  to 
the  program  could  readily  be  made  if  new  experimental  data  or  new  theoretical 
understanding  led  to  changes  in  heat  transfer  calculations  for  liquid  helium. 
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V  Conclusions  and  recommendations 

The  principal  conclusion  from  this  computer  simulation  is  that  Laketites 
show  promise  for  use  as  dielectric  coatings  at  low  temperatures.  All  of  the 
compounds  show  promise  for  improving  the  thermal  response  of  heated  super¬ 
conductors  coated  with  these  materials  relative  to  the  performance  of  super¬ 
conductors  insulated  with  GE  7031  and  other  similar  epoxies.  In  particular 
the  SC-3  and  SC-2  compounds  with  their  large  thermal  conductivities  should 
be  vigorously  developed  as  they  show  the  greatest  promise  in  keeping  super¬ 
conductors  exposed  to  large  transient  distributed  heating,  such  as  eddy  cur¬ 
rent  heating  due  to  changing  magnetic  fields,  as  cool  as  possible.  While 
not  tested  in  this  project  it  would  seem  likely  that  the  SC-3  and  SC-2  mat¬ 
erials  would  also  substantially  improve  the  ability  of  superconductors  to 
withstand  localized  heating  such  as  that  caused  by  mechanical  motion.  Super¬ 
conducting  devices  made  with  Laketite  coated  wires  would  have  greater  relia¬ 
bility  since  they  could  withstand  larger  heat  pulses  for  a  longer  time  with¬ 
out  heating  above  the  critical  transition  temperature  and  'going  normal'. 

The  additional  thermal  stability  might  also  permit  a  reduction  of  mass  of 
the  superconducting  device  by  allowing  some  of  the  copper  and/or  bronze 
usually  used  for  this  purpose  to  be  removed.  Since  the  Laketites  are  less 
dense  than  copper,  the  replacement  of  copper  by  a  Laketite  reduces  the  mass. 

These  computer  results  clearly  showed  that  the  most  important  thermal 
property  for  a  coating  to  have  is  high  thermal  conductivity.  The  SC-2  and 
SC-3  materials  have  the  greatest  promise  among  the  Laketite  materials.  One 
can  also  conclude  from  these  results,  that  thickness  is  not  a  very  important 
parameter  for  thermal  properties.  Thicker  coatings  are  slightly  advantageous 
because  they  increase  the  total  specific  heat  and  the  total  enthalpy.  On 


the  other  hand,  increasing  thickness  decreases  the  utility  of  high  thermal 
conductivity.  The  variation  of  the  superconductor  temperature  with  coating 
thickness  is  so  slow  that  the  thickness  should  be  set  by  mechanical  properties 
and  not  by  thermal  properties. 

The  most  immediate  need  for  research  work  is  for  the  experimental 
confirmation  of  the  results  of  this  computer  program.  Computer  simulations 
should  not  be  regarded  as  authorities  until  validated  with  experimental 
evidence.  Since  the  computer  program  in  this  report  assumes  that  the  outer 
edge  of  the  superconductor  is  copper,  experimental  tests  with  coated  copper 
would  be  the  best  method  of  assessing  the  validity  of  the  results  contained 
in  this  report.  The  parameters  for  specific  heat,  thermal  conductivity, 

C 

density  etc,  of  the  Laketites  were  derived  from  the  report  of  the  manufacturers 
of  the  Laketites.  This  data  should  be  verified  by  experiments  run  by  people 
not  directly  connected  with  the  manufacturer  of  these  products  to  avoid  any 
appearance  of  inflated  specifications. 

Because  of  the  promising  results  of  this  computer  simulation,  I  recommend 
a  vigorous  research  and  development  project  for  the  Laketite  materials.  Since 
the  thermal  characteristics  of  the  Laketites  are  based  on  their  ferro-electric 
nature,  the  thermal  properties  of  these  materials  should  be  tested  in  strong 
electric  fields.  The  possibility  exists  that  the  specific  heat  and  thermal 
properties  might  be  repressed  by  electric  fields.  Because  superconductors 
are  used  in  electrical  apparatus,  the  presence  of  strong  electric  fields  is 
always  a  distinct  possibility.  These  electric  tests  should  include  both 
steady  electric  fields  and  pulsed  electric  fields  such  as  might  be  present 
in  combat  areas  due  to  Electromagnetic  Pulses  (EMP).  Similarly,  the  manu¬ 
facturers  claims  that  Laketites  are  dielectrics  and  can  be  used  in  insulation. 
This  claim  should  be  verified  and  quantified  by  experimental  tests  of  the 
electrical  properties  of  these  materials  including  resistivity,  dielectric 
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constant,  A.C.  loss  parameters,  and  breakdown  resistance.  Not  only  are 
electrical  parameters  important,  but  successful  applications  of  coatings 
require  knowledge  of  mechanical  properties  at  low  temperatures  such  as 
thermal  expansion,  brittleness,  strength,  bending  radius,  fracture  strength, 
adhesion  to  the  surface  that  is  coated,  etc.  This  report  has  assumed  that 
the  surface  of  the  coating  is  mechanically  perfect,  with  no  cracks  or  voids 
and  unaffected  by  magnetic  or  electric  fields  that  are  producing  the  heat. 
The  validity  of  such  assumptions  need  to  be  tested  by  actual  experimental 
tests  of  the  materials  at  low  temperatures.  In  brief  all  relevant  pro¬ 
perties  of  the  Laketites  should  be  tested  prior  to  full  scale  development. 

I  recommend  further  development  of  the  computer  programs.  Two  parti¬ 
cular  problems  need  to  be  addressed.  One,  the  process  of  cooling  or 
'recovery'  after  heating  needs  to  be  simulated.  This  will  require  major 
modifications  of  this  program,  but  is  particularly  important  for  the 
evaluation  of  pulsed  systems,  where  the  cooling  time  is  as  important  as  the 
response  to  heating.  Two,  the  program  needs  to  be  extended  to  the  two 
dimensional  case  of  localized  'hot  spots'  since  this  is  a  more  common 
source  of  heating  of  superconductors  than  is  the  case  of  distributed  heating 
which  was  covered  in  this  report. 
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APPENDIX  I 


A  Typical  SCEPTRE  Program 


MODEL  DESCRIPTION 

MODEL  LIQUID  HELIUM  (TEMP)  (A-G) 

A  flat  surface  pool  boiling  in  saturated  liquid  helium  at  4K  (.804  ATM). 

The  program  contains  peak  nucleate  boiling  and  minimum  film  boiling 
information  specific  to  the  surface  in  use.  More  details  and  instructions 
about  the  program  can  be  found  in  the  Fortran  function  FJ  listed  in 
Appendix  II  of  the  report.  The  units  are:  V(=TEMP  IN  K) ,A(=J/S ) ,SEC, 
F(=J/K) ,OHMS(=K/WATT) 

ELEMENTS 
CL ,A-G=1. E-6 

JS,A-B=FJ ( VCL ,VCBL ,PA,PTHP ,PQP ,PTHM,PQM,1 . ) 

CBL,B-G=l.E-3 

JBL,B-C=FJ (VCL  ,VCBL,PA,PTHP ,PQP ,PTHM,PQM,2. ) 

E0,C-G=4.0 

R0,C-G=1.E4 

This  section  sets  up  the  basic  model  electric  circuit  for  the  helium, 
note  that  the  Fortran  function  FJ  given  in  Appendix  II  is  crucial. 

DEFINED  PARAMETERS 
PA=1 . 

PTHP=. 5 
PQP=. 6 
PTHM=6. 8 
PQM= . 2 

PTHP  and  PQP  describe  peak  nucleate  boiling,  respectively  they  are 
temperature  and  heat  flux,  and  PTHM  and  PQM  similarly  describe  minimum 
boiling  point.  These  parameters  are  adjustable. 

PTHL=X  T(VCL-4. 0) 

PEBL=X  E(CL*PTHL+CBL*(VCBL-4. 0) ) 

PTR=X  B(PEBL*PTHL**3. ) 

PLTR=FL(PTR,0. ) 

INITIAL  CONDITIONS 

VCL=4.0 

VCBL=4.0 

CIRCUIT  DESCRIPTION 

This  describes  the  wire  and  coating  with  W  indicating  wire  and  C  indicating 
coating.  PTW  is  the  thickness  of  the  wire  in  centimeters.  PTC  is  the 
thickness  of  the  coating  in  centimeters.  Both  parameters  are  adjustable 
as  is  PKAPC  which  sets  the  strength  of  the  kapitza  resistance. 

ELEMENTS 

J1 ,1-2=Q  J(T  JTOT(TIME) ,2. ) 

CW1 ,2-l=Q  CW(PTW,8.9,VCW1,PA,2. ) 

RW1 ,2-3=Q  R(PTW,T  KCU(VCWl) ,PA,1. ,2. ) 

02,1-3=0  J (T  JT0T (TIME ) ,2. ) 

CW2,3-1=Q  CW(PTW,8. 9,VCW2,PA,2. ) 

RW2,3-21=Q  R(PTW,T  KCU(VCW2),PA,1. ,2. ) 

CC1 ,21- 1 =Q  C(PTC,1.,T  CPC3A(VCC1 ) ,PA,3. ) 

RC1, 21-22=Q  R(PTC,T  KC3A(VCC1 ) ,PA,1.  ,3. ) 
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CC2,22-1=Q  C(PTC,1.,T  CPC3A(VVC2) ,PA,3. ) 
RC2,22-23=Q  R(PTC,T  KC3A(VCC2) ,PA,1. ,3. ) 
CC3,23-1=Q  C(PTC,1.,T  CPC3A(VCC3) ,PA,3. ) 
RC3,23-88=Q  R(PTC,T  KC3A(VCC3) ,PA,1. ,3. ) 
RK,88-99=Q  RK(PKAPC,PTS,PDT) 

H,99-1=M0DEL  LIQUID  HELIUM 
DEFINED  PARAMETERS 


PA=1 

PTW=! 0125 

PTC=.005 

PKAPC=15. 

PTS=X1(1.*VCC3) 

PTHW=Q  TH(VCWl) 

PTHS=Q  TH(PTS) 
PDT=X2(ABS(PTS-VCLH) ) 
PLSEC=FL(TIME,0. ) 

PLTHS=FL(PTS ,4. 0) 
PLEBL=FL(PEBLH,0. ) 

OUTPUTS 

IRK(QS ) ,PLOT ( PTHS(THS ) ) 

JSH(QINBL) ,PLOT(PLSEC(LOGSEC) ) 
JBLH(QOIJTBL)  ,PLOT(PLSEC(LOGSEC) ) 
PTHW(TWIRE) 

PTHS(THS) 

PTHLH(THL) 

PTHS(THS) ,PLOT(PLSEC(LOGSEC)) 
PEBLH(EINBL) 

JSH(QINBL) ,PLOT(PTHLH(THL) ) 
PLEBL(LOGEBL) ,PLOT(PLSEC(LOGSEC) ) 
PLTHS(LOGTHS ) ,PLOT(PLSEC(LOGSEC) ) 
PLTRH(LOGTR),PLOT(PLSEC(LOGSEC)) 
PLSEC(LOGSEC) 

INITIAL  CONDITIONS 
VCW1=4. 0 
VCW2=4.0 
VCC1=4.0 
VCC2=4.0 


VCC3=4.0 

FUNCTIONS 

EQUATION  J(I,N)=(I/N) 

EQUATION  R(D,K,A,GFCTR,NE)=(D*GFCTR/(K*A*NE)) 

EQUATION  C ( 0 , RHO , CP , A , NE ) = ( D*RHO*CP*A/NE ) 

EQUATION  CW(D,RHO,T,A,NE)=( (D*RHO*A* 

( (1 . 08E-5*T )+30. 6*(T/3  44.5)**3. ) )/NE) 

EQUATION  TH ( TEMP )= (TEMP-4.0) 

EQUATION  RK(C,T,DT)=(C/(T**3.*(1.+1.5*(DT/T) 

+  (DT/T)**2.  +  .25*(DT/T)**3. )))  ^  u  u  . 

JTOT  sets  the  heating  level  and  the  shape  of  the  heating  P^'se. 
is  piecewise  linear  between  discrete  points.  The  tables  of  Cp 
are  likewise  piecewise  linear  between  discrete  points. 

TABLE  JTOT 
0,0,0, .05, .001, .05 


It 

and  K 
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TABLE  KCU 

4,2.5,5,3,8,5,10,6,12,7,15,8.7,18,10,20,10.8,25,11.8, 
30,12,35,11.2,40,10,50,7.7,60,6.3,70,5.5,80,4.9,100,4.5 
TABLE  CPC3A 

4..  0018. 4. 5.. 0027. 5.. 0039. 6.. 0076. 7.. 0137. 8.. 023. 9.. 036, 

10..  055. 12.. 107. 15.. 216. 20.. 448. 30. 1.35. 40. 3. 2 
TABLE  KC3A 

4. 1.2. 5. 1.4. 6. 1.3. 7. 1.8..  8. 9.. 6. 10.. 5. 15.. 2, 

20. .  1 .25. . 1 
RUN  CONTROLS 
STOP  TIME=. 1 

INTEGRATION  R0UTINE= IMPLICIT 
TERMINATE  IF(PTHS.GT.50. ) 

The  reruns  allow  the  program  to  be  performed  again  with  a  single 
variable  changed. 

RERUN  DESCRIPTION 
FUNCTIONS 
TABLE  JTOT 
0,0,0,.l,.001,.l 
RERUN  DESCRIPTION 
FUNCTIONS 
TABLE  JTOT 
0,0,0, .2, .001, .2 
RERUN  DESCRIPTION 
FUNCTIONS 
TABLE  JTOT 
0,0,0, . 5, . 001 , .  5 
RERUN  DESCRIPTION 
FUNCTIONS 
TABLE  JTOT 
0,0, 0,1, .001,1 
RERUN  DESCRIPTION 
FUNCTIONS 
TABLE  JTOT 
0,0,0,2,-001,2 
RERUN  DESCRIPTION 
FUNCTIONS 
TABLE  JTOT 
0,0,0,3,.001,3 
RERUN  DESCRIPTION 
FUNCTIONS 
TABLE  JTOT 
0,0,0,4,.001,4 


OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOr'JOOOOOOOOOOOOO 
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FUNCTION  FJ (TL ,TBL , AREA, THLPNBS .QPNBSS 
1 ,THLMFBS,QMFBSS,JNUM) 


A  DEPENDENT  CURRENT  SOURCE  MODEL  FOR  HEAT  TRANSFER 
BETWEEN  A  SURFACE  AND  A  LIQUID  HELIUM  BOUNDARY  LAYER 
TAKES  INTO  ACCOUNT  THE  SURFACE  TEMP,  CONDITION  OF 
THE  BOUNDARY  LAYER,  VARIOUS  MODES  OF  HEAT  TRANSFER, 
TIME  HISTORY,  AND  THE  KAPITZA  RESISTANCE. 


DEFINITION  OF  SOME  TERMS  (THE  VARIABLE  NAMES  AS 
THEY  APPEAR  IN  THE  SCEPTRE  PROGRAM  ARE  ENCLOSED 
IN  <  >  BRACKETS): 


JNUM 

_ 

A  NUMBER  DESIGNATING  THE  DEPENDENT  CURRENT 

TL 

SOURCE  USING  THE  SUBPROGRAM  -  MAKES  A 

DIFFERENCE  IN  THE  HEAT  FLUX  <JS  FOR  JNUM=1. 

AND  JBL  FOR  JNUM=2.> 

TEMPERATURE  OF  THE  LIQUID  ADJACENT  TO  THE 

TBL 

- 

SURFACE.  <VCL> 

EFFECTIVE  AVERAGE  TEMPERATURE  OF  THE  BOUNDARY 

CS 

LAYER.  EQUALS  THE  VOLTAGE  ON  CAPACITOR  CBL 

(USE  TO  CALCULATE  ENERGY  IN  BOUNDARY  LAYER).  <VCBL> 

A  SMALL  CAPACITANCE  IN  THE  SURFACE  DIRECTLY 

CL 

ADJACENT  TO  THE  LIQUID.  USED  TO  GIVE  THE 

TEMPERATURE  DIRECTLY  AT  THE  SURFACE. 

A  SMALL  CAPACITANCE  IN  THE  LIQUID  DIRECTLY 

CBL 

ADJACENT  TO  THE  SURFACE.  USED  TO  GIVE  THE 

TEMPERATURE  OF  THE  LIQUID  AT  THIS  POINT.  <VCL> 

HEAT  CAPACITY  OF  THE  BOUNDARY  LAYER.  DEFINED 

AREA 

SUCH  THAT  EBLPUA=EPNBSS  WHEN  THE=THPNBSS  IN 

STEADY  STATE.  <CBL> 

SURFACE  AREA  OF  THE  OBJECT  (OFTEN  TAKEN  AS 

DLTAMAX 

— 

UNIT  AREA).  <FA> 

THICKNESS  OF  THE  FULLY  DEVELOPED  THERMAL 

THPNBSS 

BOUNDARY  LAYER  FOR  CONDUCTION  AND  CONVECTION 

HEAT  TRANSFER. 

STEADY  STATE  THETA-SURFACE  AT  THE  PNB  POINT, 

QPNBSS 

I.E.  THE  WALL  TEMPERATURE  MINUS  THE  BATH  TEMPERATURE  AT  PEAK 
NUCLEATE  BOILING  HEAT  FLUX.  <PTHP> 

STEADY  STATE  HEAT  FLUX  AT  THE  PNB  POINT.  <PQP> 

THMFBSS 

= 

STEADY  STATE  THETA-SURFACE  AT  THE  MINIMUM  FILM  BOILING 

POINT.  <PTHM> 

oooooooooooooooooooooooooooooooooooooooooc->oooooooooo 
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QMFBSS  =  STEADY  STATE  HEAT  FLUX  AT  THE  MINIMUM  FILM  BOILING 

POINT.  <PQM> 

EBNB  =  THE  VALUE  OF  DELTA  *  THL  AT  WHICH  NUCLEATE 

BOILING  BEGINS  -  DIMENSIONS  OF  CM*K  -  THIS 
NUMBER  IS  FROM  NBS  DATA  BUT  IS  CLOSE  TO  THE 
VALUE  PREDICTED  BY  THEORY  (3E-5).  <PBNB> 

TR1  =  THE  VALUE  OF  (EBLPUA*THL**3)  AT  WHICH 

TRANSITION  BOILING  BEGINS  -  THE  THL 

TERM  SEEMS  TO  COME  FROM  THE  TEMPERATURE  DEPENDENCES 

OF  THE  DENSITY,  HEAT  CAPACITY,  AND  SURFACE 

TENSION  IN  THE  HELIUM.  THE  NUMBER  FOR  TR1 

WAS  DETERMINED  FROM  NBS  DATA. 

TR2  =  THE  VALUE  OF  (EBLPUA*THL**3)  AT  WHICH  FILM 

BOILING  BEGINS  -  CAME  FROM  LOOKING  AT  THE  TEMPERATURE 
DEPENDENCE  OF  THE  MFB  POINT  AND  TRYING  TO  FIND  AN 
UNDERLYING  PRINCIPLE  WITHIN  THE  GAS  FILM 
STABILITY  ARGUMENTS  FOR  THE  MINIMUM  FILM 
BOILING  POINT.  THE  NUMBER  FOR  TR2  WAS  GUESSED 
AT,  USING  NBS  DATA  AS  A  BASE. 

EPNBSS  =  THE  ENERGY  IN  THE  BOUNDARY  LAYER  (CAPACITOR 

CBL)  AT  WHICH  PNB  OCCURS  IN  STEADY  STATE. 

EMFBSS  =  THE  ENERGY  IN  THE  BOUNDARY  LAYER  AT  WHICH  MFB  OCCURS 

IN  STEADY  STATE  HEAT  TRANSFER. 


TO  CHANGE  THIS  PROGRAM 


CHANGE  THE  DATA  POINTS  FROM  THE  QDOT/A 
VS.  THETA  SURFACE  STEADY  STATE  HEAT  TRANSFER  GRAPH 
FOR  HEAT  TRANSFER  FROM  A  PARTICULAR  SURFACE  TO  LIQUID 
HELIUM  TO  MAKE  THIS  DEPENDENT  CURRENT  SOURCE  MODEL 
REFLECT  THE  PROPERTIES  OF  THE  SURFACE. 

>NOTE:  THL  REPRESENTS  THE  LIQUID  TEMPERATURE  AT  THE 
SURFACE,  WHICH  IS  DIFFERENT  FROM  THE  SURFACE  TEMPERATURE 
DUE  TO  THE  KAPITZA  TEMPERATURE  DROP  AT  THE  SURFACE.  DO 
NOT  JUST  PLUG  IN  THE  NUMBERS  OFF  AN  ORDINARY 
Q  VS.  THS  PLOT. 

PNB:  Q<PQP>,  THL<PTHP> 

MFB:  Q<PQM>,  THL<PTHM> 

THE  HEAT  TRANSFER  EQUATIONS  USED  HERE 
SHOULD  ALSO  BE  SCRUTINIZED  FOR  APPLICABILITY 
AND  THE  DATA  CHECKED  FOR  ACCURACY. 


THE  KAPITZA  RESISTANCE  OF  THE  SURFACE  IN  CONTACT 
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C  WITH  LIQUID  HELIUM  I  MUST  BE  CHANGED  IF  THE  SURFACE 
C  CHANGES.  THE  KAPITZA  RESISTANCE  IS  ASSUMED  TO  BE 
C  OF  THE  FOLLOWING  FORM: 

C  PKAPC 

C  RK  _ _ 

C  (T**3. )*(1.-1. 5*(DT/T)+(DT/T)**2.-. 25*(DT/T)**3. ) 

C 

C  WHERE  TI  IS  THE  SURFACE  TEMPERATURE  AND  DT  IS  THE  KAPITZA 
C  TEMPERATURE  DROP  BETWEEN  TS  AND  TL.  THE  SURFACE  IS  THUS 

C  REPRESENTED  BY  ONLY  THE  CONSTANT  <PKAPC>  IN  THE  MAIN 

C  PROGRAM.  PKAPC  =91.52  FOR  GRAPHITE 
C  AND  PKAPC=12.  FOR  COPPER  SURFACES. 

C  >NOTE  THAT  A  T**(-3)  RELATIONSHIP  IS  APPROACHED  FOR 

C  SMALL  VALUES  OF  DT.  THIS  EXPRESSION  USES  THE  FULL 

C  EXPANSION  AND  NOT  JUST  THE  FIRST  TERM. 

C  REFERENCE:  NBS  TECHNICAL  NOTE  385  BY  N.S.  SNYDER 

C 

C 

C  THE  HEAT  FLUX  INPUT  TO  THE  SURFACE  IS  CHANGED  BY 

C  ALTERING  THE  TABLE  JTOT  IN  THE  MAIN  PROGRAM.  TIME 

C  AND  WATT/CM**2  COORDINATES  ALTERNATE  IN  THIS  TABLE, 

C  TIME  BEING  THE  INDEPENDENT  VARIABLE. 

C 

C 

C  IF  THE  CIRCUIT  DESCRIPTION  IS  CHANGED,  BE  SURE  TO 
C  RENAME  ALL  REFERENCES  TO  THE  SURFACE  CAPACITANCE 

C  WHICH  IS  USED  TO  PROVIDE  THE  SURFACE  TEMPERATURE. 

C  LOOK  IN:  <PTS>  AND  IN  THE  LUMPED  THERMAL 
C  RESISTANCES  NEAR  THE  SURFACE  ELEMENT. 

C  ALSO  CHECK  THE  INITIAL  CONDITIONS  SECTION. 

C 

C  A  CHANGE  IN  THE  LHE  BATH  TEMPERATURE  WOULD  PROBABLY  REQUIRE 

C  CHANGES  IN  <PTHP>,<PQP>,<PTHM>,<PQM>,  BNB,  BTB, 

C  BFB,  RHOL,  CPL,  KL,  AND  HFG.  BE  VERY  CAREFUL  ABOUT 

C  DOING  THIS. 

C 

C 

C 

C 

C  "SSHT"  STANDS  FOR  STEADY  STATE  HEAT  TRANSFER 

C  "PNBSS"  STANDS  FOR  PEAK  NUCLEATE  BOILING  POINT  IN  SSHT. 

C  "MFBSS"  STANDS  FOR  MINIMUM  FILM  BOILING  POINT  IN  SSHT. 

C  "Q"  AS  A  PREFIX  STANDS  FOR  A  HEAT  FLUX  RELATED  TERM; 

C  WHEN  FOLLOWED  BY  PNB  OR  MFB  THE  NUMBER  IS  PART  OF  THE 
C  COORDINATE  FOR  EITHER  THE  PNB  POINT  OR  THE  MFB  POINT 
C  ON  A  QDOT/A  VS.  THETA  SURFACE  SSHT  GRAPH. 

C  "TH"  AS  A  PREFIX  STANDS  FOR  A  THETA  VALUE,  THAT  IS,  THE 

C  DIFFERENCE  BETWEEN  THE  SURFACE  TEMPERATURE  AND  THE 
C  LIQUID  HELIUM  BATH  TEMPERATURE. 

C  "L"  STANDS  FOR  THE  LIQUID. 
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C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 


c 

c 

c 


c 


c 


c 


c 


c 

c 


"BL"  STANDS  FOR  BOUNDARY  LAYER  AND  USUALLY  REFERS 
TO  THE  THERMAL  BOUNDARY  LAYER  IN  THE  LIQUID  HELIUM. 

"BNB"  STANDS  FOR  BEGINNING  OF  NUCLEATE  BOILING,  AS 
EXTRACTED  FROM  EQUATION  1.13  OF  HSU  &  GRAHAM'S  BOOK 
ON  TWO-PHASE  FLUID  FLOW  AND  BOILING  HEAT  TRANSFER— 

I  EXTRACTED  THL*DELTA  AS  BEING  THE  IMPORTANT 
PARAMETER  AND  CALLED  THIS  BNB.  A  DIRECT  THEORETICAL 
CALCULATION  OF  BNB  FOR  THIS  SYSTEM  GIVES  A  VALUE  OF 
1.45E-6  CM*K.  ASSUMING  A  LINEAR  TEMPERATURE  PROFILE  AND 
CONSTANT  PROPERTIES,  THE  ENERGY  IN  THE  BL  PER  UNIT 
AREA  WOULD  BE: 

E/A  =  RH0L*CPL*(1. 45E-6 )/2  =  4.5E-7  J/CM**2 
OTHER  LIQUID  SUCH  AS  WATER  REQUIRE  HIGHER  VALUES 
THAN  THE  CRITICAL  VALUE  BEFORE  BOILING  IS  OBSERVED 
SO  I'M  NOT  SURPRISED  THAT  A  VALUE  OF  2.E-6  SEEMS 
TO  GIVE  THE  BEST  RESULTS  (4.4  TIMES  LARGER  THAN  THE 
THEORETICAL  CRITICAL  VALUE). 


PRELIMINARY  DATA  INTRODUCTION,  MODIFICATION,  AND 
CHECKING  OF  THE  INPUT  TEMPERATURE. 

REAL  KL.JNUM 

DATA  RHOL/ . 129/ ,CPL/4 . 19/ , KL/1 . 95E-4/ 

1 , DLTAMAX/2 . E-4/ ,CBL/1 . E-3/ ,CL/1 . E-7/ 

2, EBNB/2.E-6/,TRl/l.E-4/ 

3, HFG/21.47/ ,RHOV/. 01356/ 

TBATH=4.0 


FJ=0. 


PATHTR2=0 

PATHNB=0. 

PATHTB=0. 

PATHFB=0. 


EXPQNB=2. 
EXPQFB=1 . 2 


THL=TL-TBATH 

THE=TBL-TBATH 


IF(THL.LE.O. )G0T0  900 

IF( (THE.LE.O. ) .AND. (JNUM.GT. 1.5) )GOTO  900 

I F ( THE . LT . 0 . )THE“0. 


TRPNB=5. *TRI 

EPNBSS=TRPNB?(THLPNBS**3. ) 

EPNBSS  IS  4. E-3  FOR  THLPNBS=.5  AND  TRI=1 . E-4 
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TH=THLMFBS 
PATHTR2=2 
GOTO  530 
5  PATHTR2=0. 

EMFBSS=EFBSS 

TRMFB=EMFBSS*(THLMFBS**3. ) 

TR2=TRMFB*5 

EMFBSS  IS  9.57E-3  FOR  THLMFBS=6.8  AND  TR2=15.04 


HEAT  TRANSFER  MODE  DETERMINATION 

1-CALCULATE  THE  ENERGY  IN  THE  BL  PER  UNIT  SURFACE 
AREA  AND  DETERMINE  IF  FILM  OR  TRANSITION  BOILING 
EXISTS  AT  THE  SURFACE. 


EBLPUA= (THE*CBL+THL*CL )/AREA 
(EBLPUA=ENERGY  IN  THE  BOUNDARY  LAYER 
PER  UNIT  SURFACE  AREA) 

IF (EBLPUA.LE.O. )G0T0  900 

WITH  NO  ENERGY  IN  THE  BL  BUT  A  POSITIVE  THS 
HAVE  THE  CASE  OF  SIMPLE  CONDUCTION  LIMITED 
ONLY  BY  THE  KAPITZA  RESISTANCE  -  NO  CONDUCTION 
DEEP  INTO  THE  BL  IN  THIS  CASE. 

FIRST,  CHECK  TO  SEE  IF  FILM  OR  TRANSITION 
BOILING  IS  OCCURRING.  THE  PRODUCT  OF  THE 
ENERGY  IN  THE  BOUNDARY  LAYER  TIMES  THE  PEAK 
TEMPERATURE  OF  THE  LIQUID  SEEMS  TO  BE  A 
CONSTANT  FOR  THE  TWO  BOUNDARIES  OF  THE 
TRANSITION  BOILING  REGION. 

IF(JNUM.GT.1.5)GOTO  50 
EXTL=EBLPUA* (THL**3 . ) 

IF(EXTL.GE.TR2)G0T0  500 
IF(EXTL.GT.TR1)G0T0  400 

CHECK  FOR  NUCLEATE  BOILING.  EBNB  IS  AN  ENERGY 
PER  UNIT  SURFACE  AREA  (J/CM**2)  WHICH 
CORRESPONDS  TO  DELTA  *  THL  *  RHOL  *  CPL 
THE  ACTUAL  VALUE  USED  FOR  EBNB  IS  ABOUT  2 
TIMES  GREATER  THAN  THE  THEORETICAL  VALUE  AND 
WAS  TAKEN  FROM  NBS  DATA  USING  CARBON  SURFACES. 
IF ( EBLPUA . LE . EBNB )GOTO  100 
GOTO  300 

50  IF(EBLPUA.LE.EBNB)GOTO  100 
IF { EBLPUA. LE.EPNBSS)G0T0  350 
IF (EBLPUA.LT. EMFBSS )G0T0  450 


oooooooooo 
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GOTO  500 
C 
C 
C 
C 

C  MODE  I:  CONDUCTION  HEAT  TRANSFER 
C 

C  CALCULATION  OF  THE  BOUNDARY  LAYER  TEMPERATURE  PROFILE 
C  AND  THE  CONDUCTION  HEAT  TRANSFER. 

C 

C  1-TO  FIND  THE  HEAT  FLUX  OUT  OF  THE  SURFACE, 

C  ASSUME  A  TEMPERATURE  PROFILE  OF  THE  FORM: 

C  TH  =  THL*(1.-X/DLTAMAX)**EXPN 
C  AND  THEN  USE  A  RELATION  EQUATING  THE  ENERGY  IN 
C  THE  BL  TO  RHOL*CPL  TIMES  THE  INTEGRAL  OF  TEMP(TH) 

C  OVER  X=0  TO  DLTAMAX  TO  SOLVE  FOR  EXPN. 

C 

C  2- IN  THE  EVENT  THAT  TOO  MUCH  ENERGY  IS  IN  THE  BL  FOR 

C  THIS  SORT  OF  ANALYSIS,  GET  SOME  SORT  OF  ANSWER  BY 
C  GOING  TO  MODE  II:  CONVECTION  HEAT  TRANSFER 
C 

C  3-IF  SOME  CONDUCTION  OR  CONVECTION  IS  OCCURRING, 

C  FIND  THE  SLOPE  AT  . 9*DLTAMAX  AND  USE  THIS  VALUE  TO 
C  CALCULATE  A  CONDUCTION  HEAT  TRANSFER  USING: 

C  QCOND=  - KL*DTDX ( X= . 9*DLTAMAX ) 

C  WHERE  DTDX  IS  THE  SLOPE  OF  THE  TEMPERATURE  PROFILE 
C  IN  THE  BL  AT  X=.9*DLTAMAX. 

C 

100  THAVE=EBLPUA/ (DLTAMAX*RHOL*CPL ) 

THRATIO=THL/THAVE 
IF(THRATI0.LE.1.1)G0T0  200 
EXPN=THRATI0-1 . 

IF(JNUM.GT. 1. 5)G0T0  150 
X=0. 

GOTO  180 
C 

150  IF(EXPN.GT.4.)G0T0  900 
X=.9 
C 

180  DTDX=- ( EXPN*THL/DLTAMAX )* ( 1 . -X )** ( EXPN- 1 . ) 
QCOND=-KL*DTDX 
IF(PATHNB.GT.1.)G0T0  330 
FJ-QCOND*AREA 
GOTO  900 


MODE  II:  CONDUCTION  HEAT  TRANSFER  WITH  CONVECTION 
-ASSUME  CONVECTION  SETS  A  LIMIT  OFDLTAMAX  ON 
DELTA  AND  ASSUME  A  LINEAR  TEMPERATURE  PROFILE  BETWEEN 
THL  AND  0.  AT  DLTAMAX.  CALCULATE  THE  HEAT 
TRANSFER  AS  FOR  ORDINARY  CONDUCTION. 
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200  OCV=THL*KL*AREA/DLTAMAX 
IF(PATHNB.GT.1.)G0T0  330 
FJ=QCV*AREA 
GOTO  900 
C 

C  MODE  III:  NUCLEATE  BOILING  HEAT  TRANSFER 
C 

C  PEAK  NUCLEATE  BOILING  DATA  FOR  A  SPECIFIC 
C  MATERIAL  (REPRESENTED  BY  QPNBSS  AND  THPNBSS)  ARE 
C  TAKEN  INTO  ACCOUNT  WHEN  CALCULATING  THE  HEAT 
C  TRANSFER  (WHICH  GOES  ROUGHLY  AS  THE  SQUARE  OF 
C  THE  TEMPERATURE  DIFFERENCE  THS) . 

C  TRANSIENT  BOILING  HEAT  TRANSFER  MAY  BE  HIGHER 
C  OR  LOWER  THAN  STEADY  STATE  VALUES,  DEPENDING  ON 

C  PREVIOUS  HISTORY.  THIS  IS  TAKEN  INTO  ACCOUNT 

C  BY  MULTIPLIERS  (AN  ENERGY  RATIO). 

C 

C  FTR1  IS  THE  RATIO  OF  THLTR1  AT  TRANSITION 

C  INCIPIENCE  TO  THLPNBS  AT  PEAK  NUCLEATE 

C  BOILING  (WHICH  IS  NOT  PURE  NB). 

300  FTR1=.9 

THLTR1=FTR1*THLPNBS 
CQNB=QPNBSS/ (THLTR1**EXPQNB) 

C 

QNBSS=CQNB*(THL**EXPQNB) 

ENBSS=EPNBSS* (THL/THLPNBS ) 

C 

C 

C  AN  ASSUMPTION  IS  MADE  THAT  THE  HEAT 

C  TRANSFER  IN  NUCLEATE  BOILING  IS  DETERMINED 

C  BY  THE  TEMPERATURE  OF  THE  LIQUID  AT 

C  THE  SURFACE  AND  IS  AFFECTED  BY  THE  AMOUNT 

C  OF  ENERGY  ACTUALLY  IN  THE  BOUNDARY  LAYER. 

C 

IF(THL.GE.THLTR1)G0T0  310 
EXPERNB= ( 1 . - (THL/THLTR1 ) ) 
QNB=QNBSS*(EBLPUA/ENBSS)**EXPERNB 
GOTO  320 
310  QNB=QNBSS 
C 

C  IF  THE  DELTA  *  THL  PRODUCT  IS  NOT  MUCH 

C  GREATER  (4X)  THAN  THE  THRESHOLD  VALUE  FOR 

C  NUCLEATE  BOILING,  THEN  CALCULATE  HEAT  TRANSFER 

C  FOR  A  COMBINATION  OF  CONDUCTION  AND  NUCLEATE 

C  BOILING.  THE  IDEA  IS  TO  EASE  INTO  NB  AS  MORE 

C  NUCLEATION  SITES  BECOME  ACTIVE. 

320  TAU=EBLPUA/EBNB 

IF (TAU.GE.4. )GOTO  380 
PATHNB=2. 
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GOTO  100 
330  PATHNB=0. 

QN8=QNB+ (QCOND-QNB ) * (4 . -TAU  )/3. 

GOTO  380 
C 
C 

C  NUCLEATE  BOILING  HEAT  TRANSFER  OUT  OF  THE  BOUNDARY  LAYER. 

C  THE  ENERGY  IN  THE  BL  IS  USED  TO  DETERMINE  THE 

C  HEAT  TRANSFER  OUT  OF  THE  BL  BY  CALCULATING  AN 

C  EQUIVALENT  STEADY  STATE  TEMPERATURE  DIFFERENCE  (THLSS) . 

C  THIS  EFFECTIVE  TEMPERATURE  IS  THEN  USED  IN  THE  STANDARD 

C  STEADY  STATE  NUCLEATE  BOILING  CORRELATION  DEVELOPED  ABOVE 

C 

350  QNB=QPNBSS*( (  EBLPUA/EPNBSS)**EXPQNB) 

C 

380  IF(PATHTB.GT. 1. )GOTO  410 
FJ=QNB*AREA 
GOTO  900 
C 
C 
C 
C 

C  MODE  IV:  TRANSITION  BOILING 
C 

C  TRANSITION  BOILING  IS  ASSUMED  TO  BE  A 

C  COMBINATION  OF  NUCLEATE  AND  FILM  BOILING 

C  WITH  THE  PROPORTION  OF  EITHER  HEAT 

C  TRANSFER  MODE  BEING  DETERMINED  BY  EXTL 

C  AS  RELATED  TO  TR1  AND  TR2. 

C  EQUATION  (6)  FROM  KALININ,  BERLIN, 

C  KOSTYUK,  AND  NOSOVA  IN:  ADVANCES  IN 

C  CRYOGENIC  ENGINEERING,  VOL. 21  (1976), 

C  P.  273-277  WAS  USED  WITH  A  SMALL  ADDITIONAL 

C  TERM  TO  MATCH  THE  BOUNDARY  CONDITIONS  OF 

C  FRACTNB=1  AT  PNB  AND  FRACTNB=0  AT  MFB. 

C  THE  1/4  POWER  OF  EXTL  AND  THE  TRANSITION 

C  CONDITIONS  ARE  USED  TO  PRODUCE  A  RELATION 

C  WHICH  IS  LINEAR  IN  TEMPERATURE.  A 

C  DIMENSIONLESS  VARIABLE  DTSTAR  IS  USED  IN 

C  THE  EQUATION  FOR  FRACTNB  WHICH  SHOULD  BE 

C  SIMILAR  TO  THE  DIMENSIONLESS  TEMPERATURE  (5)  USED 

C  IN  THE  ABOVE  REFERENCE. 

C 

400  PATHTB=2. 

GOTO  300 
410  GOTO  500 
420  PATHTB=0. 

C 

DTSTAR= ( ( EXTL** . 25 )- (TR1** . 25 ) )/ 

1((TR2**.25)-(TR1**.25)) 

FRACTNB-O. -DTSTAR  )**10 
FRACTFB-(1. -FRACTNB) 
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QTB=QNB*FRACTNB+QFB*FRACTFB 
FJ=QTB*AREA 
GOTO  900 
C 
C 

C  TRANSITION  BOILING  -  OUT  OF  THE  BOUNDARY  LAYER 
C 

C  OTHSS  (THE  STEADY  STATE  HEAT  FLUX  IN 

C  TRANSITION  BOILING)  IS  GIVEN  BY  LINEAR 

C  INTERPOLATION  BETWEEN  THE  PEAK  NUCLEATE 

C  BOILING  AND  MINIMUM  FILM  BOILING  POINTS. 

C 

450  EXPQTSS= (ALOG(QMFBSS)-ALOG(QPNBSS) )/ 

1 (ALOG ( EMGBSS ) -ALOG ( EPNBSS ) ) 
CQTSS=QPNBSS/EPN8SS**EXPQTSS 
QTBSS=CQTSS*EBLPUA**EXPQTSS 
FJ=QTBSS*AREA 
GOTO  900 
C 
C 
C 
C 

C  MODE  V:  FILM  BOILING 
C 

500  TH=THL 

C  QFBSS  (THE  HEAT  FLUX  FOR  STEADY  STATE  FILM  BOILING 

C  AT  A  GIVEN  THS)  IS  AN  EXPONENTIAL 

C  FUNCTION  OF  THS.  THE  EQUATION  IS  AN 

C  EMPIRICAL  FIT  TO  DATA  FOR  THE  NBS  GRAPHITE 

C  SURFACES  WHICH  DID  NOT  EXHIBIT  A  MINIMUM 

C  FILM  BOILING  HEAT  FLUX  BELOW  THAT  OF  PNB. 

C  FOR  RADII  OF  CURVATURE>.5  CM. 

NPASS=0 

C 

530  QFBSS= ( . 02*TH**1 . 2 ) 

NPASS=NPASS+1 

C 

C  CALCULATE  THE  STEADY  STATE  VAPOR  LAYER 

C  THICKNESS  (DELTASS)  BY  USING  AN  EQUATION 

C  IN  YUKIO  IWASA'S  CRYOGENICS  PAPER 

C  (CRYOGENICS,  MAY  1978,  P  .273  EQ.16). 

C  UNITS  HAVE  BEEN  CHANGED  TO  GIVE  DELTA 

C  IN  CM.  THS  WAS  REPLACED  BY  THL  AND  A 

C  TERM  STATED  AS  (THS-.5)  WAS  REDUCED  TO  THL. 

C  HFG  IS  THE  LATENT  HEAT  OF  VAPORIZATION 

C  PER  UNIT  VOLUME  OF  SAT.  VAPOR  (J/CM**3). 

HFGVOL=HFG*RHOV 
DELTASS=((5. E-4*TH)/HFGV0L)+ 

1(5. 3E-5*(TH**3 . ) )/ ( 3 . *HFGVOL ) 

C 

C  THE  ENERGY  IN  THE  FILM  LAYER  IS  CALCULATED 

C  BY  USING  THE  FOLLOWING  RELATION  FOR  THE 


oooooooooo  ooooooooo  o  o  ooooo  ooooo 
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VOLUMETRIC  ENTHALPY  (J/CM**3): 

HVOL  =  .609  *  (T**(-.25)) 

WHICH  HAS  BEEN  INTEGRATED  FROM  THE  SURFACE 
OUT  TO  DELTA,  ASSUMING  A  LINEAR  TEMPERATURE 
PROFILE  BETWEEN  TL  AND  TBATH. 

EFILM=( . 812*DELTASS/TH)*( ( (TH+TBATH)**. 75)- (TBATH**.75) ) 

ENBSSFB  IS  CALCULATED  ASSUMING  THAT 
THE  SAME  ENERGY  IS  ADDED  AS  FOR  REGULAR 
NUCLEATE  BOILING  WITH  A  HEAT  FLUX 
EQUAL  TO  QFBSS. 

ENBSSFB=EPNBSS*( (QFBSS/QPNBSS)**(1 ./EXPQNB) ) 

EFBSS=EFILM+ENBSSFB 
IF(PATHTR2.GT. 1 . )G0T0  5 
ERAT 10= EBLPUA/ EFBSS 

IF(JNUM.GT. 1. 5 )G0T0  550 

ASSUME  A  LINEAR  TEMPERATURE  PROFILE 
SO  THAT  THE  TEMPERATURE  SLOPE  AT  THE  SURFACE  IS 
DIRECTLY  PROPORTIONAL  TO  THE  RATIO  OF 
(EFBSS/EBLPUA) 

SINCE  THS  IS  GIVEN.  THIS  ALSO  ASSUMES 
THAT  ANY  TEMPERATURE  DEPENDENCE  OF  THE  FLUID 
PROPERTIES  WILL  DROP  OUT  IN  THE  RATIO. 

QFB=QFBSS*( EFBSS/EBLPUA) 

IF(PATHTB.GT. 1 . )GOTO  420 

FJ=QFB*AREA 

GOTO  900 


FILM  BOILING  -  OUT  OF  THE  BOUNDARY  LAYER 

CALCULATE  AN  EQUIVALENT  THLSS.  ASSUME 
THAT  FILM  BOILING  HEAT  TRANSFER  OUT 
OF  THE  BOUNDARY  LAYER  IS  DETERMINED 
BY  EBLPUA.  THIS  IS  TO  MEET  THE 
CONDITION  THAT  QINBL=QOUTBL  WHEN 
EBLPUA=EFBSS. 

550  IF{ (ERATIO.GT. .9999) 

1 . AND . ( ERAT 10 . LT . 1 . 0001 ) )GOTO  580 
TH1=TH 
C 

TH=TH1*(EBLPUS/ EFBSS) 

I F (NPASS.LT. 10 )GOTO  530 
C 

TH=TH1*( ( EBLPUA/ EFBSS )+ 1 . )/2. 

I F ( NPASS . LT . 50 )GOTO  530 
C 

TH=TH1*( (EBLPUA/EFBSS)+2. )/3. 

IF (NPASS.LT. 100 )G0T0  530 
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FILMED 


